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ABSTRACT 



In order to study the nucleate pool-boiling performance 
of R-114 and R- 114-oil mixtures from enhanced evaporator 
tube surfaces, an experimental apparatus was designed, 
constructed and instrumented. The evaporator was made of a 
T-shaped Pyrex glass container. Boiling occurred from a 
smooth, hard-copper tube, 15.9 mm (5/8 in.) in outer diam- 
eter, 12.7 mm (1/2 in.) in inside diameter and 431.8 mm (17 
in. ) in length. The tube was heated using a cartridge 
heater, and it was instrumented with 8 thermocouples to 
measure the wall temperature. A Hewlett-Packard 3497A data 
acquisition/ control unit and a 9826A computer were used to 
collect and process data. The condenser was cooled by an 
ethylene glycol-water mixture, which was maintained at about 
-17 °C by means of an R- 12 refrigeration system. Nine data 
runs were completed to de-bug the experimental apparatus and 
to check for reproducibility. During .all data runs, espe- 
cially at higher heat fluxes (greater than 10 kW/m^), large 
temperature variations were observed along and around the 
active boiling length of the test tube. The data were 
compared with data found in the literature and reasonable 
agreement was obtained. 
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tube outside surface area of active boiling section 
cross-sectional area of tube 

tube outside surface area of non-boilng section 
conversion coefficients for thermocouple 
specific heat 
diameter 

tube inside diameter 
tube outside diameter 
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latent heat of vaporization 

current 

current reading by AC Current Sensor 
thermal conductivity of liquid 
thermal conductivity of copper 
active boiling tube length 
non-boiling length of the test tube 
Nusselt number 
pressure 
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tube outside wall perimeter 
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heat- transfer rate from boiling surface 
heat-transfer rate through one unboiling end 
heat -transfer rate from cartridge heater 
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average wall temperature at the thermocouple 
location 

calibration bath temperature 
temperature of the thermocouple location 
saturation temperature of the boiling liquid 
outer wall temperature of the boiling test tube 
voltage across the cartridge heater 
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volumetric thermal expansion coefficient 
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I. INTRODUCTION 



A. BACKGROUND 

In order to reduce the size and weight of air- 
conditioning systems and water chillers, the United States 
Navy recently has considered two different approaches: 1) to 
change the type of the refrigerant from R-11 to R-114, and 
2) to use enhanced evaporator tube surfaces, such as Gewa-T, 
Thermoexcel-E , and High-Flux. Schematic drawings and photo- 
micrographs of these enhanced surfaces at various magnifica- 
tions are given in Figures 1.1 , 1.2 and 1.3. A combination 
of R-114 and enhanced surfaces may lead to significant 
improvements in system efficiency as well as to smaller and 
lighter components in future Navy applications. 

B. ADVANTAGES OF USING R-114 

The latent heat per unit volume of R-114 is 50 BTU/ft^ 
at 70 °F. As a result, the heat carried by each cubic foot 
of R-114 is nearly twice as much as that carried by R-11 
[Ref. 1] at the same temperature. In other words, the 
volume of R-114 vapor circulated per ton of refrigeration is 
about half that for R-11. Inherently, the size of the system 
components can be reduced significantly using R-114 rather 
than R-11. R-12 and R-22 in fact can transfer the heat more 
efficiently than R-114, but they require higher pressure 
ranges. However, operating pressures for R-114 are moderate, 
similar to R-11, such that light-weight components may be 
desirable and less-experienced personnel may be involved in 
the installation and maintenance. While the refrigerants in 
general are nonflammable, R-114 and R-12 are rated to have 
the lowest toxicity among all refrigerants [Ref. 2]. Based 
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on the above discussion, R-114 is probably the most suitable 
refrigerant for use on board submarines and ships. 

Further, R-114 is considerably more stable toward heat 
than R-12. All standard ferrous and nonferrous metals, 
except zinc and magnesium, and aluminum alloys containing 
appreciable amounts of zinc and magnesium, may be used with 
R-114 [Ref. 3]. The physical properties of R-114 and its 
comparison with the other common refrigerants are summarized 
in Appendix A, pressure- enthalpy diagram of R-114 is 
provided in Appendix B and the range of application of 
refrigerants are provided in Appendix C [Ref. 4]. 

C. ENHANCED BOILING TUBE SURFACES 

The evalution of enhanced surface geometries that 
promote high-performance nucleate boiling has been surveyed 
by Webb [Ref. 5]. Two basic types of enhanced boiling 
surfaces are employed on the outer surface of tubes; porous 
coatings and reentrant grooves. The High-Flux surface is a 
sintered, porous metallic matrix bonded to the base surface. 
The Gewa-T and Thermoexcel-E surfaces have reentrant - type 
nucleation sites formed by cold working an integral finned 
tube wall. See Figures 1.1, 1.2 and 1.3 for details. 

The key to the high performance of these surfaces is 
attributed to three factors [Ref. 5] : "(1) a pore or reent- 
rant cavity within a critical size range, (2) interconnected 
cavities, and (3) nucleation sites of a reentrant shape. 
When the cavities are interconnected , one active cavity can 
activate adjacent cavities. It appears that the dominant 
fraction of the vaporization occurs at a very thin liquid 
film within the subsurface structure. The reentrant cavity 
shape provides a stable vapor trap, which will remain active 
at very low liquid superheat values.” 
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The pool boiling heat - 1 rans f er performance of R-113 from 
the above-mentioned commercially-available enhanced surfaces 
has been investigated by Marto and Lepere [Ref. 6]. They 
reported that for R-113, the enhanced surfaces showed a two 
to tenfold increase in the heat- transfer coefficient when 
compared to a plain tube. The High-Flux surface was most 
effective over a broad range of heat fluxes, whereas the 
Gewa-T surface was only effective at high heat fluxes (near 
100 kW/m^). The Thermoexcel-E surface showed similar gains 
in the heat- transfer coefficient to that of the High-Flux 
surface below 10 kW/m^ . They also reported that the degree 
of superheat required to activate the enhanced surfaces was 
less than the plain tube, and was sensitive to initial 
surface conditioning and fluid properties. 

D. NUCLEATE POOL BOILING PERFORMANCE OF PURE R-114 

The heat - 1 rans fer performance of pure R-114 in .the 
evaporator section of an air-conditioning unit can be esti- 
mated using one of the suggested commonly-jused correlations 
for nucleate pool boiling, which are listed by 
Chongrungreong and Sauer [Ref. 7]. However, because of the 
lack of complete understanding of the nucleate boiling 
mechanism, these correlations should be compared with exper- 
imental data if it is possible. So far, the literature 
reveals a limited amount of nucleate boiling experimental 
investigations using R-114. Happel [Ref. 8], using a hori- 
zontal, pure-nickel tube with a surface roughness of 0.43 
micrometer, investigated the boiling characteris tics of pure 
R-114 at various pressures. The effect of pressure on 
nucleate boiling heat transfer in R-114 was also investi- 
gated experimentally by Nishikawa and his co-workers 
[Ref. 9]. They operated over a reduced pressure (P/P^ ) 
range from 0.03 to 0.98, and demonstrated that the heat 
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transfer coefficient of R-114 increases remarkably at high 
pressure near the critical point. Their results showed that 
the effect of surface roughness changes with pressure and 
the effect gradually vanishes as pressure approaches the 
critical point. Using a horizontal smooth copper tube, 
Henrici [Ref. 10] investigated boiling heat- transfer 
performance of R-114 over a heat-flux range from 0.1 to 100 
kW/m^ and over a pressure range from 37 to 252 kN/m^ . 

Stephan and Mitrovic [Ref. 11] investigated the perform- 
ance of pure R-114 from a Gewa-T enhanced surface. Boiling 
occurred outside a horizontal bundle of these tubes. They 
reported that the position of the tube within the bundle 
does not effect the heat - 1 rans f er coefficient considerably. 
The reason for this effect was explained in the following 
manner: ’’The flow pattern around a tube seems to be of very 
weak influence and the heat - transfer process is mainly 
governed by the bubble formation and the two-phase flow 
inside the T-shaped channel. The top part of the fins, which 
were rolled to form the T, practically did not directly 
contribute to the heat - transfer process because, even at 
maximum heat fluxes , vapour bubbles were not produced 
there . ” 

Marto and Hernandez [Ref. 12] also investigated nucleate 
pool-boiling characteris t ics of a Gewa-T surface in R-113. 
Their experimental results showed a disagreement with 
Stephan and Mitrovic 's above explanation. They showed that 
’’the liquid-vapor motion may not always be confined to flow 
circumferentially within the channels formed by the T-shaped 
fins as previously postulated. Surface imperfections on the 
tips of the Gewa-T fins may cause bubble nucleation from 
these sites. The presence of these active sites may further 
enhance the performance of the Gewa-T surface.” 

The factors which appear to affect boiling heat - transfer 
performance of a single tube are the physical properties of 
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refrigerant, boiling pressure, tube diameter, surface condi- 
tion of the tube, and maybe the "submergence effect" 
(resulting from the hydrostatic liquid head) on the evapora- 
tion process. Experimental observations have shown that 
changes in magnitude of these properties and conditions 
significantly affect pool-boiling heat- transfer ’ performance . 
[Ref. 7]. 

E. NUCLEATE BOILING PERFORMANCE OF REFRIGERANT- OIL MIXTURES 

Most of the refrigeration systems work with oil- 
lubricated, hermetically sealed compressors, which are 
cooled by the working fluid. Therefore, the fluids being 
circulated within these refrigeration systems generally are 
not pure refrigerants, but are refrigerant-oil mixtures. 
The saturation temperature of a ref rigerant - oil mixture is 
higher than that for the pure refrigerant at the same pres- 
sure. Also, both viscosity and surface tension increase 
when oil is added to refrigerants. Other thermodynamic prop- 
erties of refrigerants, such as density and specific heat, 
are also significantly changed by the addition of oil. Thus, 
adequate knowledge on the influence of oil on the boiling of 
refrigerants is very important for the successful design of 
refrigeration systems. The literature reveals only very 
little information on the boiling performance of R- 114-oil 
mixtures [Ref. 10]. In particular, heat- transfer perform- 
ance of the enhanced evaporator surfaces in ref rigerant - oi 1 
mixtures is a poorly-understood subject area. Since these 
advanced boiling surfaces contain a high density of pores, 
the presence of oil can generally result in a greater 
performance degradation than on more- convent ional , smooth 
surfaces . 

Various studies which are surveyed and listed by Jensen 
and Jackman [Ref. 13] have shown that, in general, nucleate 
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boiling heat - t rans f er coefficients decrease with increasing 
oil concentration in R-11, R- 12 , R-22, R-113 and R-114. For 
oil concentrations less than 3%, however, Dougherty and 
Sauer [Ref. 14] demonstrated that the heat-tansfer coeffi- 
cient may be increased slightly for R-ll-oil and R-113-oil 
mixtures. They attributed this trend to an increase of 
foaming of the mixtures. This result also has been • confirmed 
by Sauer, Gibson and Chongrungreong [Ref. 15] using R-12-oil 
mixtures, and by Henrici [Ref. 10] using R-ll4-oil mixtures. 
However, in a recent experimental study using R-113-oil 
mixtures, Jensen and Jackman [Ref. 13] did not observe such 
an enhancement with low oil concentrations. The model of 
the governing processes involved in boiling of refrigerant- 
oil mixtures has been given by Stephan [Ref. 16] as follows: 
"Obviously the reduction in heat transfer in a mixture is 
due to the fact that the liquid near the vapour- liquid 
interface of growing bubbles becomes denuded in the more 
volatile components. Thus, a concentration gradient exists 
between the interface and the bulk of the liquid, and the 
more volatile components have to diffuse through the liquid 
in order to reach the interface. A mass- transfer resistance 
thus is introduced, which is not present when boiling single 
component s . " 

The effect of the oil viscosity on heat transfer in 
boiling of ref rigerant - oi 1 mixtures has been considered by 
Sauer, Gibson, and Chongrungreong [Ref. 15]. Oil viscosity 
is believed to significantly contribute to the changes in 
boiling characteristics of ref rigerant -oil mixtures. The 
heat- transfer coefficient is greatest for the highest 
viscosity of oil for the same temperature difference. 

The effect of oil on boiling of R- 12 from a Gewa-T tube 
was investigated by Stephan and Mitrovic [Ref. 17]. They 
reported that the presence of oil in R- 12 reduced the heat- 
transfer coefficient almost over the entire range of heat 
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fluxes and mass fractions. They observed an increase in 
heat- transfer coefficient with less than 5% mass fraction of 
oil only at high heat fluxes (such as higher than 20 kW/m^). 

The effects of oil on the heat- transfer performance of 
refrigerants in a tube bundle have been summarized by Arai , 
et al. [Ref. 18] in the following manner; "The presence of 
oil in a shell- and- tube evaporator is believed to have two 
opposing effects on the heat- transfer coefficient of the 
evaporator tube bundles. Firstly, many people suspect that 
the oil left out from the evaporating refrigerant tends to 
choke the surface structure, thereby lowering the heat- 
transfer coefficient. Secondly, from the experience of 
conventional heat exchangers, it has long been recognized 
that the presence of oil enhances foaming on the upper part 
of the tube bundles. This foaming increases the vigor of 
refrigerant movement around the tubes, so that heat transfer 
is enhanced." 

F. THESIS OBJECTIVES 

Based on the foregoing discussion, it is clear that 
little information exists in the literature to enable the 
successful design of R-114 evaporators, especially when 
using enhanced boiling surfaces, in the presence of lubri- 
cating oils. The objectives of this research effort were 
therefore to: 

1. Design, construct and instrument an experimental 
apparatus in order to investigate the heat - trans f er 
performance of R-114 and R-114-oil mixtures from 
enhanced evaporator tube surfaces , 

2. Successfully operate the system to produce repeatable 
data, and 

3. Obtain smooth- tube baseline data for R-114 to be used 
as a standard for the future data on advanced 
surfaces . 
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Figure 1.1 Schematic and Optical Hicrograph of 
Section of the Gewa-T Surface 12-OX). 



a Cross 
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Tunnel 



Figure 1.2 



Schematic and Optical Micrograph of the 
Thermoexcel - E Surface (40X) . 
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Figure 1 . 3 



Schematic and 
High Flux 



Scanning Electron Micrograph of the 
Surface ( 500X ) . 
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II. DESCRIPTION OF EXPERIMENTAL APPARATUS 



A. OVERVIEW OF THE SYSTEM 

An overall schematic representation of the experimental 
apparatus is shown in Figure 2.1, and a photograph is shown 
in Figure 2.2. The apparatus consists of a Pyrex-glass 
evaporator with the boiling test tube 1, a Pyrex-glass 
condenser 2, an R-114 liquid reservoir 3, an oil reservoir 
4, a graduated oil cylinder 5, a vacuum pump 6, a cooling 
section consisting of a 1/2-Ton R-12 refrigeration unit, a 
tank filled with water- ethylene glycol mixture 7, and a pump 
8. The evaporator 1, condenser 2, oil-reservoir 4, 
graduated-oil cylinder 5 and R-114 reservoir 3, are placed 
in a Plexiglas chamber 9 , where the temperature can be main- 
tained at lower values (up to 10 °C) relative to ambient 
temperature. This chamber also served as a safety barrier in 
case of an inadvertent overpressurization of the glass 
apparatus . 

The vapor rising from the evaporator is condensed in the 
condenser 2 and the liquid is then fed back into the evapo- 
rator through valve V5 . In this manner, the specific 
enthalpy of condensation ("latent heat") of R-114 vapor is 
removed in the condenser. The evaporator of the R-12 
refrigeration system is made of a copper coil, which is 
immersed in the water-ethylene glycol tank in order to main- 
tain the temperature of the mixture at about -17 °C. An 
8-GPM, turbine-type pump 8 transports the water- ethylene 
glycol mixture through the control valve VC to the 
condenser. 

The oil can be added into the evaporator 1 using the 
graduated cylinder 5 by opening valves VI and V4 . The 
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graduated cylinder is supplied with oil by reservoir 4. The 
valves V5 , V6 and V7 are provided for emptying the evapo- 
rator prior to dismantling it, for example, for changing the 
evaporator tube. Also, fresh R-114 can be added to the 
reservoir while valve V5 is closed and V6, V7 are open and 
by providing R-114 vapor through V8 and VIO . 

The test tube was placed horizontally in the evaporator 
after a 240 VAC electrical cartridge heater of 203.2 mm (8 
in.) in length had been inserted in its middle section. The 
whole system was connected to the vacuum pump through V8 and 
VIO in order to remove noncondensible gases. The relief 
valve was set to a gage pressure limit of 138 kN/m^ (20 
psi), which is 50% less than the manufacturer- recommended 
working pressure limit of the Pyrex glass tees. All connec- 
tions were assembled with Swagelok fittings sealed with 
Teflon ferrules to ensure leak tightness. Also, all valves 
used in the system, except control valve VC and by-pass 
valve V9, were made of stainless steel. Copper tubes of 
9.5-mm (3/8 in.) diameter were used for all piping, except 
in the oil section, where 6.35-mm (1/4 in.) tubes were used. 

B. BOILING TEST SECTION 

1 . Evaporator 

The evaporator is a T-shaped container made of 
Corning Pyrex glass. Figure 2.3 shows a schematic of the 
evaporator with dimensions. The selection of Pyrex glass 
offers several advantages: it is corrosion- res is tant and 
transparent and has a smooth interior surface (this mini- 
mized nucleate boiling at the inner surface of the 
container) compared to metals, and is more stable with 
temperature and pressure variations compared to ordinary 
glass. The operating pressure limit of the evaporator is 
guaranteed by the manufacturer up to a gage pressure of 200 
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kN/m^ (30 psi). Each end of the Pyrex-glass evaporator was 
fitted with a cast-iron flange and a gasket. A detailed 
sketch of the cast-iron flange is shown in Figure 2.3. Two 
aluminum flanges, 210 mm in diameter and 12.7 mm in thick- 
ness, were bolted to the cast-iron flanges. Thus, all 
fittings were connected through the aluminum flanges to the 
Pyrex-glass evaporator. 

2 . Test Section 

A schematic drawing of the test section is shown in 
Figure 2.4. The boiling test tube was held in place by two 
Teflon bushings , which were attached to the aluminum flanges 
at both ends of the T-shaped evaporator. Four studs were 
used in order to attach the Teflon bushing on the aluminum 
flange. The Teflon bushing and both ends of the test tube 
were sealed by means of two rubber 0-rings in order to 
achieve leak tightness. An additional 0-ring was placed 
between the aluminum flange and the Teflon plug at each end. 

3 . Boiling Tube 

A smooth, hard-copper tube, 15.9 mm (5/8 in.) in 
outer diameter, 12.7 mm (1/2 in.) in inside diameter and 
431.8 mm (17 in.) in length, was used to provide baseline 
data for the nucleate pool boiling heat - transfer coefficient 
of pure R-114. The heater was a 1000-Watt 240-Volt 
stainless - steel cartridge, 6.35 mm (1/4 in.) in outer diam- 
eter and 203.2 mm (8 in.) in length. The cartridge heater 
was inserted into a copper sleeve, which was 6.35 mm (1/4 
in.) in inside diameter, 12.7 mm (1/2 in.) in outside diam- 
eter and 203.2 mm (8 in. ) in length. In order to provide a 
uniform heat flux, the cartridge heater and the copper 
sleeve were soldered together. Pulido [Ref. 19] showed that 
(using a similar type of heater- element construction) the 
heater element provides a uniform heat flux along the 
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circumferential direction. The heater element used during 
this thesis was assumed to have similar characteristics. 
The cartridge heater and the copper sleeve were then 
inserted as a unit into the middle portion of the test tube 
using a tight mechanical fit. (Note: A shrink-fit process 
is essential to minimize the thermal contact resistance. 
However, owing to the time constraint in this thesis, the 
shrink-fit process was delayed and a slide fit with a clear- 
ance of less than 0.01 mm was used. Thus, the resulting tube 
was used mainly for the purpose of obtaining a successfully- 
operating apparatus). The active boiling length of the test 
tube was 203.2 mm (8 in.) in the middle portion of the 
boiling tube. Tn order to compute the actual average heat 
flux in the heated portion of the tube, a suitable correc- 
tion was applied for the natural convection created at both 
ends . See Heat Flux Calculation section in Chapter 4 for 
details . 

To measure the wall surface temperature of the 
boiling tube, 8 thermocouples were inserted into 8 grooves 
which were machined on the outside of the copper sleeve. As 
shown in Figure 2.5, these thermocouples were located at 
different axial and circumferential locations. Four thermo- 
couples (1,2, 5, 6) were placed with 90-degree separation 
around the copper sleeve to measure the temperature distri- 
bution at the mid cross section of the active length of 
boiling tube. The thermocouple's wire was glued down at 
several places along the thermocouple's channel using Epoxy. 
The longitudinal temperature distribution on the active 
boiling section was also measured by the symmetrically- 
located 4 thermocouples (3, 4, 7, 8). The exact locations of 
the thermocouples and dimensions of the thermocouple grooves 
are given in Figure 2.5. All thermocouple grooves were 
axially machined from the location of the thermocouple hot 
junctions to the nearest end of the copper sleeve. 
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C. CONDENSER SECTION 



The condenser was also a T-shaped container made of 
Corning Pyrex glass. It was identical to the evaporator. The 
position of the condenser can be seen in Figure 2.1. R-114 
was condensed on a vertical copper coil, which was inserted 
in the Pyrex-glass condenser. The copper coil was fabricated 
into a 76.2 mm (3 in.) diameter coil using 9.5 mm (3/8 in.) 
copper tube. The active condensation length was estimated to 
be 4.5 m (15 ft). 

The top portion of the condenser was connected to a 
portable, mechanical vacuum pump to remove noncondensible 
gases from the apparatus . The bottom of the condenser was 
also connected to the evaporator via valve V5 in order to 
return the condensed R-114 liquid to the evaporator. The 
cooling liquid, i.e., water-ethylene glycol mixture, entered 
the top portion of the condenser, through the copper coil 
and left the condenser from the bottom to return to the 
water- ethylene glycol tank. The condenser was placed verti- 
cally and connected to the vapor outlet of the evaporator 
using L-shaped aliminum tube, 50.8 mm (2 in.) in diameter. 
The maximum vapor velocity of R-114 vapor through the 
aluminum tube was found to be about 0.5 m/sec. A Bourdon 
gage with a range of absolute vacuum to a gage pressure of 
1030 kN/m^ (150 psi) and a relief valve which was set to 138 
kN/m^ (20 psi) were placed on the L-shaped aluminum tube. 

D. OIL ADDING SECTION 

To study the boiling performance of R-114-oil mixtures, 
a cylindrical aluminum reservoir, 152.4 mm (6. in) in diam- 
eter and 152.4 cm (6 in.) in height, and a glass oil 
cylinder were installed above the evaporator. The relative 
positions of the oil reservoir and the oil cylinder are 
shown in Figure 2.1. The oil cylinder was 355 mm in length 
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and had a diameter of 25.4 mm. This cylinder was specially 
ordered to achieve a resolution of 0.5 ml. The oil cylinder 
was connected to the oil reservoir through valves V3 and V2 . 
The addition of oil into the evaporator can be achieved 
through VI by gravity after balancing the pressure of the 
oil cylinder with that in the evaporator by opening valve 
V4. 

E. COOLING SECTION 

1. Water- Ethylene Glycol Mixture Tank 

In order to store the water- ethylene glycol mixture, 
a special tank was manufactured. The total volume of the 
tank was 0.154 m^ (0.48 m x 0.48 m x 0.66 m) and it was made 
of 12 . 7-mm- thick Plexiglas sheet. All sides of the tank were 
glued together with methylene- chloride solution. The joints 
were held together with small screws for extra strength. The 
low thermal conductivity of Plexiglas was especially suited 
to minimize heat transfer (from room to water- ethylene 
glycol mixture) through the tank walls. The tank was placed 
on the floor and all sides were insulated with 22 mm (7/8 
in.) thick rubber insulation sheets. The cooling mixture 
contained 49 liters (13 gal.) of ethylene glycol and 94 
liters of (25 gal.) distilled water. The freezing point of 
this mixture was about - 25 °C. 

2 . R- 12 Refrigeration Plant 

A 1/2-Ton R-12 refrigeration plant was installed to 
cool the water- ethylene glycol mixture. Figure 2.6 shows a 
schematic of the R-12 refrigeration plant. It consists of a 
compact-type air-cooled condenser, a compressor, a receiver, 
a filter-drier unit, a pressure regulator, a pressure- 
control switch and a thermostatic expansion valve. The 
evaporator of the R-12 refrigeration plant was constructed 
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using a 9.5 mm (3/8 in.) copper tube, which was immersed in 
the water- ethylene glycol tank. The temperature of the 
water- ethylene glycol mixture is controlled by both a ther- 
mostatic expansion valve and a pressure control switch. The 
R-12 refrigeration plant was adjusted to keep the tempera- 
ture of the cooling liquid at about -17 °C. 

3 . Pump and Cont ro 1 Valve 

An 8 GPM, 115 VAC Burks turbine- type , positive- 
displacement pump was installed on the floor and the 
suction side of the pump, 25.4 mm (1 in.) in diameter, was 
directly coupled to the water- ethylene glycol tank. Cooling 
liquid was pumped from the tank to the condenser through the 
control valve VC. Also, a by-pass valve V9 was placed before 
the control valve VC on the discharge line. The use of the 
by-pass line served two important purposes: 1, it avoided 
overloading of pump 8 in the event value VC is completely 
closed, and 2, it provided proper mixing for the "warm” 
stream returning from the condenser 2. See Figure 2.1 for 
the positions of control valve VC and by-pass valve V9 . The 
by-pass valve was adjusted and set permanently so that a 
sufficient range of flow rates can be achieved through valve 
VC to satisfy proper cooling of the condenser. 

F. R-114 RESERVOIR 

An aluminum cylindrical reservoir, 228 mm (9 in.) in 
diameter and 254 mm (10 in.) in height, was placed verti- 
cally between the evaporator and condenser in order to store 
R-114 as a liquid. The liquid level of the R-114 can be 
observed by means of a sight hose attached on the reservoir 
with proper fittings. R-114 reservoir was connected to the 
vapor line through valve V7 and to the liquid line through 
valve V6. See Figure 2.1 for arrangement of the reservoir. 
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G . CHAMBER 



An aluminum frame (1.07 mx 0.51m x 0.61 m) was 
constructed to locate all the parts of the apparatus , except 
the cooling section. All four vertical sides of the frame 
were covered with 12.7 mm (1/2 in.) thick Plexiglas sheets 
and both left and right sides were provided with hinges to 
enable easy access to the components of the apparatus. 
Aluminum and plywood plates were used to cover the bottom 
and top sides of the frame, respectively. The valve bodies 
of VI through V8 were placed inside of the front Plexiglas 
sheet with the valve stems penetrating the sheet. Thus, the 
valve handles were accessible from outside of the Plexiglas 
box. The whole frame was placed above the water- ethylene 
glycol tank with aluminum support so that the system was 
very compact. 

One of the main advantages of this chamber is that the 
temperature surrounding the evaporator can be reduced rela- 
tive to the ambient temperature. Also, in case of emergency, 
the thick Plexiglas chamber would provide a safety barrier 
to personnel and equipment. 

H. INSTRUMENTATION 

1 . Power Measurement 

A 240 Volt AC source was used as the power supply, 
and it was adjusted by a variac in the range of 0-260 Volt 
and 0-8 Ampere according to the desired heat flux at the 
surface of the boiling tube. Power input to the boiling tube 
was measured with an AC current sensor and an AC-DC true 
R.M.S converter. The AC current sensor was connected to the 
input line of the heater in series and the AC-DC true R.M.S 
converter was connected in parallel. Figure 2.7 shows a 
schematic representation of the power-measurement devices . 
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Both the AC current sensor and the AC-DC true R.M.S 
converter were connected to the data acquisition/ reduct ion 
system. 

2 . Temperature Measurement 

Various temperatures were monitored throughout the 
system to include: 

1. Boiling tube wall (8 thermocouples) 

2. Liquid temperature (one thermocouple) 

3. Vapor temperature (two thermocouples), and 

4. Water- ethylene glycol mixture temperature (one ther- 
mocouple) 

The locations of the wall thermocouples in the sleeve of the 
boiling tube are shown in Figure 2.5. The liquid and vapor 
thermocouples were inserted into the two specially- 
manufactured thermocouple wells. Figure 2.8 shows a sche- 
matic of these thermocouple wells. While the stainless- steel 
portion minimizes (owing to low thermal conductivity) errors 
resulting from the axial conduction of heat from the 
surrounding, the copper tip helps minimize the temperature 
drop from the area being measured to the thermocouple loca- 
tion (owing to the high thermal conductivity of copper) . 

All the temperature measurements were accomplished 
by 0.245 mm (30 gage) copper- cons tantan thermocouples. Each 
thermocouple measurement was read directly by a 
Hewlett-Packard 3497A data acquisition system, which was 
controlled by a Hewlett-Packard 9826 computer. Each thermo- 
couple was scanned for 0.8 seconds and twenty readings were 
averaged to obtain a more accurate measurement. 

A total of five thermocouples were calibrated. Two 
thermocouples were made from the beginning of a spool of 
copper- cons tantan wire; one was made from the mid portion 
and two were from the end. It was assumed that the proper- 
ties of the copper and constantan do not change along a 
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given section of wire for any 
were calibrated by the method 



given spool. All thermocouples 
described in Appendix D. 
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Figure 2.1 Schematic of the Boiling Test Apparatus. 
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Figure 2.2 



Photograph of Overall System. 
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a) Coming Pyrex Glass Evaporator (D x El = 402x51 mm, 
L = 17 8 mm , L]_ = 12 7 mm ) 










INSERT 

t 

b) Cast Iron Flange and Gasket (d 3 _= 190mm, do = 210mm, 
= 14mm , A = 21*^) 



Figure 2.3 Schematic of the Pyrex Glass Evaporator. 
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Figure 2.6 Schematic of R-12 Refrigeration Plant. 
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Figure 2.7 Schematic of the Power Measurement. 
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Figure 2.8 Sketch of a Thermocouple Well. 



III. DATA ACQUISITION/ REDUCTION 



A. DATA ACQUISITION AND STORAGE 

A Hewlett-Packard 3497A automatic data acquisition/ 
control unit was used to read temperatures from the thermo- 
couples and to read current and voltage values of the 
cartridge heater from the AC current sensor and the AC-DC 
true RMS converter, respectively. A Hewlett-Packard 9826A 
computer was used to control the data acquisition/ control 
unit and to analyze and store data. Figure 3.1 shows a 
photograph of the data acquisition/ contro 1 unit and 
computer . 

Information was entered through the keyboard to prompt 
the data acquisition/ control unit to automatically scan each 
channel. Channel assignments are listed in Table I. These 
raw data were immediately processed and a hard-copy printout 
was provided. Also, these data were transferred to a 
computer disk under a user- specif ied file name to keep a 
permanent record. The ability to store raw data directly 
enabled these data to be reduced at any time and allowed 
flexibility for changes to data- reduct ion software. 

B. DATA REDUCTION 

Following data acquisition for each data point, results 
were computed according to the stepwise procedure outlined 
in the next section, and then printed on a Hewlett-Packard 
2671G thermal printer. Heat flux versus wall superheat (the 
temperature of the boiling surface minus the liquid tempera- 
ture) were also stored in a user- specif ied plot data file 
for subsequent plotting using the subprogram PLOT. A 
Hewlett-Packard 7470A plotter was also interfaced with the 
computer . 
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C. STEPWISE DATA- COLLECTION AND SOLUTION PROCEDURE 

1. Select tube type (all dimensions of the boiling test 
tube are included) . 

2. Set desired heat flux and saturation temperature of 
the boiling liquid. Wait for steady-state conditions 
(See Chapter 4 section B for details). 

3. Scan all channels listed in Table I (thermocouples, 
AC current sensor and AC-DC true RliS converter read- 
ings ) . 

4. Save these raw data in a user- specif ied file. 

5. Convert these raw data readings to corresponding 
units (temperature in °C, current in Ampere, voltage 
in volt ) . 

6. Compute the heat - transfer rate from the cartridge 
heater (See Chapter 4, Section C for details). 

7. Compute the average wall temperature of the boiling 
test tube and calculate the wall superheat of this 
data set (See Chapter 4, Section C for details). 

8. Compute the physical properties of R-114 using given 
correlations at film temperature (See Appendix H for 
details ) . 

9. Compute the natural- convection heat - transfer coeffi- 
cient of R-114 from non-boiling ends of the test 
tube . 

10. Compute heat losses from non-boiling ends. 

11. Calculate the heat flux from boiling test tube to the 
boiling liquid. 

12. Calculate boiling heat- transfer coefficient of the 
R-114 from test tube. 

13. Store the heat flux versus wall superheat values for 
each data set in a user- specif ied plot file. 

14. Use the subprogram PLOT and plot the data run. 



43 



Appendix E shows a listing of computer program DRPRl and 
subprogram PLOT and Appendix F shows an example of represen- 
tative data run. 





TABLE I 


HP 3497A 


Channel Assignments 


Channel 


As s ignment 


25 - 32 


Boiling tube wall temperature 


33 


Liquid temperature 


34-35 


Vapor temperature 


36 


Temperature of cooling liquid 


62 


AC-DC true RSM converter 


63 


AC current sensor 
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Figure 3.1 Photograph of Data Acquisition/Reduct ion System. 







IV. EXPERIMENTAL PROCEDURE 



A. PREPARATION 

1. Pressure Test of the Apparatus with Ai r 

Upon assembling all the components, including the 
boiling test tube, the experimental apparatus was pressur- 
ized with air up to 100 kN/m^ (15 psi). A soap-bubble test 
was first carried out to locate any leaks from each compo- 
nent of the apparatus. All detected leaks were succesfully 
and systematically fixed at the end of each pressure test. 

2. Pres sure Test of the Apparatus with R- 114 

A second pressure test was performed using R-114 
vapor. The apparatus were charged with R-114 vapor in the 
following manner. The apparatus was evacuated to 25 in. Hg. 
by a portable mechanical vacuum pump through valves Vll and 
V8 (See Figure 2.1 for the configuration of valves). Upon 
closing valve Vll, a vapor outlet of the R-114 supply 
cylinder was connected to the apparatus by means of valve 
VIO . Thus, the pressure of the apparatus was increased up to 
the saturation pressure of R-114 at that ambient temperature 
(ambient temperature was 21 °C and the corresponding satura- 
tion pressure (gage) was 86 kN/m^ (12.5 psi)). An Automatic 
Halogen Leak Detector, TIF 5000, was used to detect R-114 
leakage from the apparatus. The sensitivity of this detector 
is 1/2 oz. per year, 3 PPM max. A few small leaks were 
observed and successfully isolated. 

3 . Charging the Apparatus with R- 114 

Upon completing the pressure tests, the evaporator 
of the apparatus was filled with R-114 liquid up to a marked 
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position (20 mm above the boiling test tube) with the 
following procedure (See also Figure 2.1). 

1. The temperature of the water-ethylene glycol tank was 
reduced to about -17 °C by the R-12 refrigeration 
unit . 

2. The apparatus was evacuated to about 29 in Hg vacuum 
.by the vacuum pump (the valves Vll and V8 were open 

and VIO was closed). 

3. Upon closing Vll, the vapor outlet of the R-114 
supply cylinder was connected to VIO. 

4. The cooling section pump 8, delivered the cooling 
liquid from the water- ethylene glycol tank to the 
condenser through the control valve VC. 

5. The R-114 vapor condensed on the condenser coils and 
the liquid was collected in the evaporator by 
gravity, while V5 was open and V6 was closed. 

6. Pump 8 was then stopped when the desired liquid level 
of the evaporator was achieved. 

7. The pressure of the apparatus was then allowed to 
increase up to the saturation pressure corresponding 
to the ambient temperature. 

B. NORMAL OPERATION 

The following procedure was established to obtain the 
heat - 1 rans f er coefficient of R-114 from the smooth copper 
test tube: 

1. The test tube was immersed in the pool (20 mm below 
the liquid level). 

2. The R-12 refrigeration unit was operated 24 hours in 
advance in order to reduce the temperature of the 
water- ethylene glycol tank to about -17 °C. 

3. The data acquisition/ control unit, computer, power 
supply and cooling section pump were switched on. 
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4. 



Any accumulated noncondensible gases were evacuated 
by the portable mechanical vacuum pump through valve 
V8 and Vll (see Figure 2.1). 

5. The iteractive computer program DRPRl was operated 
and the desired heat flux and saturation temperature 
were given to the program as reference values. 

6. The cartridge-heater voltage was then adjusted by 
variac in order to maintain the desired heat flux. 
The desired and actual heat fluxes were compared 
continiously by the computer program until they 
agreed to within 2 percent. 

7. The actual versus desired saturation temperature of 
the liquid R-114 was also monitored by computer 
program. The amount of cooling liquid, which was 
being circulated through the copper coil in the 
condenser, was regulated by control valve VC in order 
to obtain nearly constant saturation temperature 
(±0.2 °C) at a given heat flux. 

8. For all the consecutive settings during both 
increasing and decreasing heat fluxes, the boiling 
was allowed to stabilize for five minutes at each 
power setting and raw data (thermocouple readings, AC 
current sensor readings and AC-DC true RMS converter 
readings) were recorded in a user- specif ied file. 

9. At a given power setting and a saturation tempera- 

ture, the following processed data were recorded as a 
printout: wall temperatures of the boiling section, 

liquid bulk temperature, vapor temperature, net heat 
flux, temperature of the water- ethylene glycol 
mixture, wall superheat and heat - trans f er coefficient 
of the R-114. 

10. For each data set, the above procedure beginning with 
step 5 was repeated. 
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11. For each data series, seven different heat fluxes (1, 
2, 5, 10, 20, 50 and 95 kW/m^ ) were selected. 



C. HEAT -FLUX CALCULATION 



According 


to the 


description 


of 


the boiling 


tube 


construction. 


the cartri 


dge heater. 


which 


was inserted 


into 


the boiling tube is the 


heat source. 


The 


variations in 


heat 


flux are made 


through 


different vo 


It age 


settings of 


the 



variac . 



= VI (1) 



where : 



Qh = 


heat- trans f er rate 


from the cartridge 




heater (W) 




V = 


Voltage across the 


cartridge heater 




element (volts) 




I = 


current through the 


heater element (amps) 



This circuit is connected in 
sensor and in parallel with an 
Each of these two sensors produce 
of 0-10 V, which is automat 
acquis it ion/ control unit. The cal 
sensors give the actual values as 



series with an AC current 



AC-DC true RMS converter. 


s a DC output in 


the 


range 


ically read by 


the 


data 


ibration equations 


of 


these 


follows : 







V = 25Vs (2) 

I = Is (3) 
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where : 



V 



s 



voltage reading by AC-DC true RMS 
converter (volts) 



I 



s 



current reading by AC current 
sensor (volts) 



The total length of the boiling test tube was 355.6 mm, 
while the active boiling length; i.e., the length of the 
cartridge heater, was 203.2 mm. The geometry of the boiling 
test tube is shown schematically in Figure 4.1. 



Figure 4.1 Geometry of the Boiling Test Tube. 

Natural - convect ion heat transfer was assumed to occur at 
both ends of the test tube (the distance ) . This part of 
the tube was considered to be a straight fin of uniform 
cross section. For this cylindrical geometry: 



EPOXY BOILING COPPER 
PLUG TUBE SLEEVE 




A 



c 



7T (D^ - Dp/4 



(4) 
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where : 



(5) 



A 



s 




A = cross-sectional area of tube (m^ ) 
c 

D = tube outside diameter (m^ ) 

D = tube inside diameter (m^) 

A = tube outside surface area of 
s 

nonboiling section (m^ ) 
p = tube outside wall perimeter (m) 

P = ^ D Q 

L = non-boiling length of the test tube (m) 



It was assumed that the temperature at the base of the 
straight fin was equal to the average wall temperature of 
the active boiling section. The average outer wall tempera- 
ture of the active boiling section was calculated using the 
radial conduction equation from the thermocouple position to 
the surface of the tube. 



8 




n = l 



( 6 ) 



Two= Tkvg- Di)/(2 L ) ) (7) 



where : 



= temperature of the thermocouple 
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location (°C) 

average wall temperature at the 
thermocouple locations (°C) 

D-j^ = diameter at the position of the 
thermocouple (m) 

D 2 - outer diameter of the boiling tube (m) 

L = active boiling tube length (m) 

= thermal conduct ivitiy of the copper (W/m.K) 

= outer wall temperature of the boiling 
test tube (°C) 

= temperature at the base of straight fin (°C) 

The temperature distribution of a fin of uniform cross- 
sectional area is given by [Ref. 20] with the following 
as sumptions : 

1. One-dimensional conduction in the x direction. 

2. Steady-state condition. 

3. Constant thermal conductivity. 

4. Negligible radiation from the surface of the fin. 

5. No heat generation in the fin. 

6. Uniform convective heat - transfer coefficient over the 
fin surface. 

d^6 / dx^ - m^© = 0 (8) 

where : 



0(x) - T(x) - 



(9) 
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and 



= h p / k^Ac 



( 10 ) 



where : 

T ^ = saturation temperature of fluid (°C) 

S a u 

h = convective heat - transfer 
coefficient (W/m^.K) 

Now, assuming negligible heat loss from the thin tube 
tip (at X = 

(d9 / 0 (11) 

and also at the base: 

^ -k A (dT/dx) = -1<^A (d0/dx) (12) 

^ c c x=0 ^ ^ x=0 



where : 



Qb = heat - 1 rans f er rate from the base (W) 

Using the above boundary conditions (Eq. 11, 12), the 

temperature distribution along the straight fin can be 
expressed as : 

(G/G]^) = ( Cosh(m(L^- x)) / Cosh(mL^) ) (13) 
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and the fin heat - transfer rate is: 



where : 



Qf 



Qb 



•yj h P Ac 9bTanh(mLu) 



(14) 



Qp = heat - 1 rans f er rate through one non-boiling 
length of the test tube (W) 



The average difference between the wall temperature and 
the liquid temperature may now be determined from the 
following equation: 



'u 



(e/0b) = 
e = 
0 = 



( /(Cosh(m(L^- x) )/Cosh(mL^) ) dx ) / (15) 

^0 



(9b/ mLu ) Tanh(mL^^ 



■wo 



sat 



where: 

“ '^at“ average difference between outer 

wall temperature and liquid saturation 
temperature (K) 



A free convection correlation stated by Churchill and 
Chu [Ref. 20] was applied for the average Nusselt number on 
a horizontal cylinder: 



0.387 Ra 



1/6 

D 



Nu 



D 



0. 60+- 



(1 + (0.559/Pr)®''^®) ) 



(16) 



- 5 — 12 

10 < Ra^ < 10 

o 
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The average Nusselt number is: 



Nu 



D 



o 



h D 

“X 



o 



where : 



k 



thermal conductivity of R-114 (W/m.K' 



Solving for h: 



h = 



k 1 

^ : 0 . 60 + 

o 



0.387 

o 



r”'u + (0.559/Pr)"/l*'l®/2M 



Pr = - 

a 



Ra, 



g g <^Wn- ^SAT> Pq 

V a 



6 = 



L 

p 



where : 



Pr = 


Prandtl number 




V z: 


kinematic viscosity 


(in' /s ) 


a = 


thermal diffusivity 


(in'/s) 



( 17 ) 



(18) 



(19) 



( 20 ) 



( 21 ) 
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■ average Rayleigh number 

g = gravitational acceleration (m/s^) 

S = volumetric thermal expansion 
coefficient (1/K). 

P = density of R-114 liquid (kg/m^) 



Now, substitution of Equations 10, 15, 20 into Equation 

18 results in: 



h = ^ <; 0 . 60 + 0 . 387 



gSD^e Tanhir-^) ^ LuO 
^ o B { ) 

vaL,(^) 

[1 + (0 . 559/Pr) 



( 22 ) 



V 

Equation 22 is solved for h by an iterative technique 
within a range of precision of 0.001. Knowing the value of 
the natural convection heat -transfer coefficient along the 



non-boiling ends, 
from Equation 14, 



the total heat-loss rate is calculated 
and the heat- transfer rate through the 
boiling section is obtained by subtracting the total heat- 
loss rate from the total heat- transfer rate: 






(23) 



Q = Q„- 2xQ, 



(24) 
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Finally, the heat flux from the boiling surface is; 



q = Q / 



( 25 ) 



where : 



Ab = 
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V. RESULTS AND DISCUSSION 



A. OUTLINE OF THE DATA RUNS 

Using the procedure outlined in Chapter 3, nine data 
runs were completed, primarily to de-bug the experimental 
apparatus for its successful operation. For all data runs, 
the same smooth copper tube, described in Chapter 2, was 
used as the boiling surface. Each data run consisted of 
seven different heat fluxes (1000, 2000, 5000, 10000, 20000, 
50000 and 95000 W/m^) with a specified saturation tempera- 
ture (0, 10 or 20 °C). All the raw data were stored in a 
user- specif ied file, named "WHxx” , where ”W" indicates 
Wieland tube, "H" indicates hard copper and "xx” indicates 
the data run number. Also, the heat- transfer coefficient 
versus heat flux data were stored in a plot file named 
"Pxx" . A summary of these data runs is given in Table II. 

B. LONGITUDINAL AND CIRCUMFERENTIAL TEMPERATURE VARIATIONS 

During all data runs, especially at high heat fluxes 
(greater than 10 kW/m^), large temperature variations were 
observed along the active boiling length of the test tube. 
Figure 5.1 shows the measured wall temperatures at each of 
the thermocouple locations, while Figure 5.2 shows the 
temperature variations in the circumferent ial direction (at 
the axial mid point) at two different heat fluxes (2 and 92 
kW/m^ ) . Despite the uniform heat flux provided by the 
heater, the temperature measurements showed considerable 
variations in both axial and circumf erent ial directions. 
The reason for these temperature variations may be explained 
as follows: As explained in Chapter 2, during the manufac- 
turing of the boiling test tube, the copper sleeve was 
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inserted into the boiling test tube without any shrink-fit 
process. The slip-fit process (with a diameter clearance of 
about 0.01 mm) used in this study was believed to have 
caused considerable thermal contact resistance at the inter- 
face between the boiling tube and the copper sleeve. Also, 
it is quite possible that this contact resistance was not 
uniform along both the axial and circumf erential directions . 
The circumferential temperature distribution may also be 
attributed to a circumferential variation in boiling heat- 
transfer coefficient. During the high heat flux boiling 
regime, it was observed that the rate of bubble formation at 
the top of the horizontal boiling tube was lower than at the 
rest of the boiling tube; i.e., the boiling heat - t rans f er 
coefficient at top of the boiling tube was relatively lower 
than the other locations . This poorer performance at the top 
of the tube may be attributed to the fact that this portion 
was not receiving sufficient amount of fresh liquid to 
replace the liquid that evaporated. The obstruction of the 
liquid flow could be caused by the growing vapor blanket 
from the bottom to the top of the tube. Based on this obser- 
vation, it is expected that the temperature at top of the 
boiling tube could be greater than the other locations . As 
shown in Figure 5.2, location 5, which represents the top of 
the horizontal boiling tube, measured the highest tempera- 
ture difference as expected from the above-mentioned obser- 
vat ion . 

However, positions 1 and 2 showed a sizable discrepancy 
at high heat flux, and this was due most likely to the non- 
uniform contact resistance. Also, a comparison of thermocou- 
ples .2, 7 and 8 show a similar discrepancy. It is clear 
therefore, that the shrink-fit process must be used in the 
future . 
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C. PLOT ANALYSIS OF NUCLEATE BOILING REGIME 



Figure 5.3 shows a typical nucleate pool-boiling 
performance curve of the smooth copper tube in R-114. The 
observed behavior of this process is analyzed by studying 
the different regions of the boiling curve. 

From point A to point B, a continuous increase in wall 
superheat (T^q observed when heat flux is increased. 
Also, during the experimental runs, no bubbles were observed 
along the boiling tube in this region of the curve. This 
region corresponds to the natural- convection process. 

From point B to point C, a reduction in wall superheat 
is observed, while the heat flux continuously increased. 
While the middle portion of boiling tube showed bubbles, the 
rest of the boiling tube showed no bubbles during the exper- 
iment in this region. This region is known as the mixed 
boiling region, where transition from natural convection to 
nucleate pool boiling heat transfer takes place. The wall 
superheat continued to decrease until numerous nucleation 
sites became active (point C). 

After point C, the wall superheat starts to increase 
with increasing heat flux as shown by region C to D. A very 
high density of bubble formation was observed along the 
active boiling section of the test tube (only a few bubbles 
were observed on the nonboiling sections of the test tube). 
See Figure 5.4. This observation confirmed the validity of 
the natural convection model assumed for data reduction. 
However, the R-114 liquid returning to the evaporator from 
the condenser caused a mild turbulance in the entrance 
region of this return line in the evaporator. Inherently, 
this affects the natural - convect ive heat - trans f e r perform- 
ance of nonboiling ends of the test tube. 

When heat flux is gradually decreased, the curve follows 
a different path (from D to E) as shown in Figure 5.3. This 
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is due to the existing, stable nucleation sites remaining 
active for a wide range of heat fluxes. Similar analysis can 
also explain Figure 5.5, which represents the same data run 
but on a different basis (heat - trans f er coefficient versus 
heat flux) . 

D. REPRODUCIBILITY TEST OF THE APPARATUS 

In order to test the reproducibility of the experimental 
apparatus, two runs were performed on two different days 
(data run number 5 and 7 in Table II) at the same saturation 
temperature. Figure 5.6 and Figure 5.7 show the comparison 
of these two data runs. It can be seen from these figures, 
that there is very good agreement between these two 
different runs. This agreement shows the ability of the 
apparatus to reproduce data runs revealing successful 
operation . 

E. BOILING PERFORMANCE OF SMOOTH COPPER TUBE IN R-114 

In order to test the validity of the data taken from the 
present experimental apparatus (despite the presence of 
errors owing to the thermal contact resistance mentioned 
earlier) , an attempt was made to compare the present data 
with data found in the literature. For this purpose, data 
run number 5 (10 °C saturation temperature with decreasing 
heat-flux rate) was compared with two sources in the litera- 
ture as discussed below. 

1. Comparison with Chongrungreong- Sauer Correlation 

A correlation proposed by Chongrungreong and Sauer 
[Ref. 7] for the nucleate boiling performance of refrige- 
rants and ref rigerant - oil mixtures is compared with run 
number 5. This correlation is based on a dimensional anal- 
ysis and data from various sources; [Ref. 7] , [Ref. 14] , 
[Ref. 15], [Ref. 22] and [Ref. 23]. 
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( 26 ) 



X IpjlSi'i f 

lO.OIoSS. 



f,, I S79 

I /'vi 



where : 



h = heat - t rans f er coefficient (W/m^.K) 



Q = heat rate (W) 



A = surface area (m^ ) 



D = characteristic dimension of heated surface (m^ ) 



= viscosity of saturated liquid (gr/m.s) 

I 

= latent heat of vaporization (W.s/g) 



C = specific heat of liquid (J/g.K) 

Jr 

1cl= thermal conductivity of liquid (W/m.K) 

P = pressure (atmospheres) 

^ = volume fraction of pure refrigerant 
= liquid density (g/cm^) 

= vapor density (g/cm’) 

Eq ' n (26) agrees very well with theirdata [Ref. 7] 



for R-11 as well as with Stephan's data [Ref. 22] for oil- 
free refrigerants, R-11, R-12, R-13, R-21, R-22, R-113 and 
R-114 with errors smaller than 16 percent. Figure 5.8 pres- 
ents a comparison of the current experimental data with eq ' n 
(26). As can be seen from this figure, the current data are 
in excellent agreement with this predictive equation. (The 
analysis given in Appendix I shows that the uncertainty in 
the boiling heat - t rans fer coefficient for this experimental 
investigation is about 20%. ) This agreement was somewhat 
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unexpected since the present data includes some kind of 
contact resistance in it. If this contact resistance were 
known, and were removed from the data, a higher boiling 
coefficient would result. 

2. Comparison with Data of Henrici 

The nucleate boiling heat - trans f er coefficient of 
R-114 and R-114-oil mixtures from a horizontal, smooth 
copper tube (30 mm in outer diameter) has been experimen- 
tally investigated by Henrici [Ref. 10]. Since Henrici 's 
data were taken from a tube with a larger diameter than the 
current boiling tube diameter, his data were revised to 
include this difference. For this purpose, a semi- emprical 
correlation to show the effect of tube diameter, as devel- 
oped by Cornwell et al . [Ref. 21], was used: 

Nu = C Re^ ' ^ (27 ) 



where : 



Re = (q D) / (h^gP ) and Nu = (h D) / k 
C = 150 for refrigerant 



Using this relationship, Henrici 's data were revised 
for a diameter of 15.875 mm (5/8 in.) used in this 
experiment . 



h(D=15.875 mm) = h(D=30 mm) (30/15.875)'^' 
h(D=15.875 mm) = 1.24 h(D=30 mm) 
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A line representing the revised Henrici’s data is 
plotted in Figure 5.9 together with the unrevised data line. 
It can be seen that the boiling heat-transfer coefficient 
increases by about 25 percent when the diameter is decreased 
from 30 mm to 15.875 mm (50-percent reduction in diameter). 

Figure 5.10 provides Henrici's natural convection 
data and Henrici's revised boiling data in comparison with 
the present experimental data. As expected, the natural- 
convective heat - transfer coefficient is lower than the 
boiling heat- transfer coefficient. The boiling heat- transfer 
coefficient from Henrici's data is higher than the coeffi- 
cient found in this study. The reason for this may be 
explained as follows. In Henrici's experimental investiga- 
tion, the measurement of the surface temperature has been 
done using thermocouples outside of the boiling surface. 
Unfortunately, he did not describe in detail how he measured 
the wall temperatures. If the thermocouples were not 
totally imbedded in the wall, they may produce incorrect 
wall temperature values. If a portion of a thermocouple is 
in direct contact with the liquid, this thermocouple will 
measure a value lower than the actual wall temperature and 
such an error would decrease the wall superheat , thus 
resulting in a higher heat - transfer coefficient. Also, if 
the contact resistance effect is removed from the present 
data, then the heat transfer coefficient would increase 
toward the Henrici data. 

F. EFFECT OF PRESSURE 

The effect of pressure on the nucleate-boiling heat- 
transfer performance of R-114 was investigated for reduced 
pressures (V/'P^ ) of 0.0268, 0.0391 and 0.0556 (data run 
numbers 5, 8 and 9, respectively in Table II). Figure 5.11 
shows the relative increase in boiling heat - t rans f er 
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coefficient with increase of pressure as expected. In order 
to calculate the pressure effect on the boiling performance, 
a simple equation developed by [Ref. 7] was used. 



where : 



h = 6.17 (q) 



0.5 5 



(<l» 



3.6 5 



. 0,2 4 



(28) 



h = heat- transfer coefficient (W/m^.K) 
q = heat flux (W/m^ ) 

0 = volume fraction of pure refrigerant 
P = boiling pressure (atmospheres) 

Using this equation, the pressure effect on the boiling 
heat transfer coefficient can be written as follows: 



h(Pi) = hCPj) (P^ /Pj )'• 



2 4 



According to the above equation, the effect of pressure 
from 127 kN/m^ to 87 kN/m^ decreases the boiling heat- 
transfer coefficient by about ten percent, which agrees 
reasonably well with the current experimental data. 
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TABLE II 
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Figure 5.1 Temperature Variations on the Boiling Tube. 
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Figure 5.3 Typical Nucleate Pool Boiling Curve for R-114 
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Figure 5.4 A Photograph of the Boiling Tube During High Heat Flux Operation 
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Figure 5.5 Typical Nucleate Pool Boiling Heat Transfer Coefficient Curve for 
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Figure 5 . 6 Reproducibility of the Apparatus . 
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Figure 5.7 Reproducibility of the Heat-Transfer 
Coefficients . 
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Figure 5.8 A comparison of the Boiling Performance with the 
Correlation of Chongrungreong and Sauer (Ref. 7). 
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Figure 5.9 Effect of Tube Diameter on Boiling Performance. 
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Figure 5.10 Comparison of Current Data with Data of Henrici 

(Ref. 10) 
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Figure 5.11 The Effect of Pressure on Boiling Performance. 



VI. CONCLUSIONS 



1. An experimental apparatus to study the nucleate- 
boiling heat - t rans f er performance of R-114 from a 
single horizontal tube has been designed, constructed 
and instrumented. 

2. An analysis scheme to establish the nucleate boiling 
heat - transfer coefficient of R-114 has been outlined, 
and data acquisition/ reduction programs have been 
written and tested. 

3. The experimental apparatus was successfully operated 
in the saturation temperature range of 0 °C to 20 °C 
and in the heat flux range of 1 to 95 kW/m^. 

4. The test runs performed on different days on the 
apparatus showed repeatabilty within ± 1 percent , 
revealing successful operation of the apparatus. 

5. The uncertainty in wall superhet (about 20%) drasti- 
cally affected the boiling heat - trans f er coefficient. 
This can be attributed to the large temperature vari- 
ations measured along and around the boiling tube 
because of a non-uniform contact pressure between the 
copper sleeve and the boiling tube. 

6. The data taken on this experimental apparatus have 
been compared with the experimental correlation for 
refrigerant and refrigerant-oil mixtures developed by 
Chongrungreong and Sauer [Ref. 7]. The present data 
agreed to within 5 percent with this correlation. 

7. The current data were also compared with Henrici's 
[Ref. 10] experimental nucleate boiling data (from 
smooth horizontal copper tube in R-114). The present 
data lie as much as 50 percent below Henrici's data. 
This disagreement may be attributable to the thermal 
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contact resistance present in the boiling tube, or to 
the errors in Henrici's wall temperature 
measurements . 
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VII. RECOMMENDATIONS 



1. To prevent considerable temperature variations along 
and around the boiling tube, a shrink- fit process 
must be applied in order to have more uniform contact 
pressure between the copper sleeve and the boiling 
tube . 

2. To determine the effect of contact resistance 
(between copper sleeve and boiling tube) on the 
nucleate-boiling heat- transfer performance of a hori- 
zontal tube, a series of experimental investigations 
should be performed with different contact pressures. 

3. To maintain constant electrical power input into the 
boiling tube, a voltage regulator should be added to 
the input line. 

4. To keep the temperature of the water- ethylene glycol 
mixture at a constant value, the 1/2- ton capacity 
R-12 refrigeration unit should be replac-ed with a 
2 - ton uni t . 

5. The positive-displacement pump should be replaced 
with a cent ri fugal - type pump, in order to deliver 
lower flow rates without reaching an electrical over- 
load limit . 

6. To check the saturation temperature of the boiling 
liquid corresponding to the measured saturation pres- 
sure, a more accurate pressure transducer should be 
used . 

7. The active length of the cartridge heater should be 
carefully measured or obtained from the manufacturer. 



80 



APPENDIX A 



PHYSICAL PROPERTIES 



OF FREON FLUOROCARBON 



COMPOUNDS 







"FREON** 11 


"FREON" 12 


"FREON** 113 


"FREON" 114 


Chemical Tormuit 




CCl.F 


CCI F, 


CCLF CCIF, 


CCIF, CCIF, 


Molecular Weight 


137 37 


120 92 


187 38 


170 93 


Boiling Point at 1 atm 


•c 


23 82 


- 29 79 


47 57 


3 77 




•F 


74 87 


-21 62 


117 63 


38 78 


Freerinf Point 


•c 


-111 


- 158 


-35 


-94 




•F 


- 168 


-252 


-31 


- 137 


Critical Temp^raturt 


•c 


198 0 


1120 


214 1 


145 7 




•F 


388 4 


233 6 


4174 


294 3 


Critical Pressurt 


atm 


43 5 


40 6 


33 7 


.32 2 




Ibs/sq In abs 


639 5 


596 9 


495 


473 2 


Critical trolumt 


cc/moJ 


247 


217 


325 


293 




ca ft/lb 


0 0289 


0 0287 


0 0278 


0 0275 


Critical Density 


f/cc 


0 5M 


0 558 


0.576 


0 582 




Ibs/cu ft 


34 6 


34 8 


360 


36 32 


Density, Liquid 


g/cc 


1 476 


1 311 


1 565 


1 456 


at 25'C (77*0 


fbs^cu ft 


92 14 


81 84 


97 69 


90 91 


at 54.5-C (130T) 


g/cc 


1 40 


1 19 


1 49 


1 16 




Ibs/cu ft 


87 40 


74 28 


93 01 


84 90 


Density, Sit'd Vapor 


g/l 


5 86 


6 33 


7 38 


7 83 


at Boiling Point 


fbs/cu ft 


0.367 


0 395 


0 461 


0 489 


Specific Heat, liquid 












(H^at Capacity) 


cal/(g) (*C) 


0 208 


0 232 


0 218 


0 243 


at 25'C (77*F) 


Blu/(lb) (*F) 










Specific Heat. Vapor. 












at Const Pressure (I atm) 


cal/(t) CO 


ni49« 38-C 

(100‘F) 








at 25'C (77*F) 


Btu/(ib) CF) 


0 145 


0 170 


Specilic Heat Ratio ol Vapor 












at 2S*C and 1 atm 


Cp/Cy 


1 137 A 38*C 

1.13/ 9 (100"F) 


1.137 


. 080 


1 084 


Heat of Vaporisation 


cal/g 


43 10 


39 47 


35 07 


32 51 


at Boiling Point 


Btu/ib 


77.51 


71.04 


63 12 


58 53 


Thermal Conductivity at 25*C (77*F) 










Blu/(hf) (ft) fF) 
Liquid* 




0 0506 


0 0405 


0 0434 


0 0372 


Vapor (1 atm) 




0.00451 


0.00576 


00014 


00060 


(Data Irom ASHRAC In most cues) 






(0 5 atm) 




Viscosity at 25'C (77'F) 












Liquid 


centipoisa 


0 430 


0 214 


0.68 


0 36 


Vapor (1 atm) 


centipoisa 


0 0107 


00123 


OHIO 


00112 


(Dala Irom ASHRAE In most cates) 






(0 1 aim) 




Surface Tension at 25*C (77*F) 


18 


9 


173 


12 > 


dynei/cm 












Refractiva lnde.i of Liquid 




1.374 


1.287 


1 354 


1 288 


at 25*C (77*F) 












Relative Dielectric Strength ol Vapor at I atm 


3.71 


2.46 


3 9 (0 44 aim) 


3 34 


and 25*C (77*F) (nitrogen 


*= 1) 










Dielectric Constant 












Liquid 




2 28 ® 29*C 


1 003^ 


2 41 ® 25'C 


2 26 O 25''C 


Vapor (1 atm)*« 




1.0036 9 24-C^ 




1.0043 <9 26 8*C 


SoluMllty of “Freon*’ In Witar 








0017 




at 1 atm and 25*C (77’F) 


wt % 


0.11 


0 028 


(Sai n Pres) 


0013 


Soluhility of Water In "Freon" 












at 25*C (77*F) 


Wt % 


0.011 


0 009 


0 011 


0 009 


Solubility Parameter ( ) 




75 


61 


72 


62 


Kauri Butanol Vaiua (KB) 




60 


18 


32 


12 


Toxicity* 
Threshold Limit 


ppm (y/») 


1000 


lOOO 


1000 


1000 


Vaiua (TLV) 


mg/m’ 


5600 


4950 


7600 


7000 
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PRESSURE- ENTHALPY DIAGRAM OF R-114 



§1118 I 9 883S98 
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APPEN DIX C 

APPLICATIONS OF "FREON" FLUOROCARBON COMPOUNDS 



The tabfe t>efow ii intended to provide a general view ol the range of applications and is not all inclusive. 
For specialized applications or more detail, please make specilic request. 



F tuorocarbon 


Ref rigerants 


A«rotol Propeltantslbl 


Solvents. Blowing Agents. 
Fire E xtingoishinii. 
Dielectric F luidt end Other Uies 


"Freon" 14 


— 


— 


— 


j "Frenn" 23 


Component of "Freon" 503 
azeotropel^l 


— 


— 


"Freon" 13 


Specialty low temperature 
applications, usually in "cascade" 
systems 


— 


— 


"F reon" 1 1 6 


— 




Dielectric fluid 


"Freon" 138 1 


Intermediate between "Freon" 13 
and "Freon" 22 lor medium to low 
temperature applications Not 
extensively used. 


— 


Efficient fire exiinguishani 
(Halon* 1301) especially suited lor 
automatic protection ol materials 
subject to water damage and of areas 
occupied by personnel 


"Freon" 22 


Household and commercial 
relnperation and air conditioning 
applications. Permits use ol smaller 
equipment. Component ol 
azeotfopes<<*l. 


— 


— 


"Freon" 1 15 


Used as an azeotrope component in 
"Freon" 502*«l . 


Accepted as a food propellant by the 
FDA. this material is well suited for 
food aerosols and (mds use in fat 
emulsion food whips Good foam 
stability with absence ol odor or taste. 


Dielectric lluid. an economic 
replacement for "Freon" 116 in 
most dielectric applications. 


"Freon" 12 


Most widely used relrigerant in 
household, automntive and 
commercial relngeration and air 
condiiioning sysienis Also as a 
component of azeotropesf''f and. 
rn high pnnty lorm ("Freon" 
freezant) approved as a direct 
contact Ireezing agent for foods. 


Most widely used high pressure 
proiiellant fur non food use Blends 
with "Freon" 1 1 and "Freon" 1 14 
are widely used. 


Blowing agent lor foamed plastics 
applicationslfii. Dielectric gas 


"Freon" 1 14 


In large industrial process cooling 
and air conditioning systems using 
multr stage ceninlugal compresson. 


Low pressure propellant, alternative 
to "freon" 1 1. having poorer 
solubility properties and less odor. 
Especially used m personal products 


Blowing agent for foamed plastics. 


"Freon" 21 




— 


Heat transfer fluid 


"Freon" 1 1 


Widely used in centrifugal 
compressors for industrial and 
commercial atr conditioning systems 
and for industrial cooling of process 
water or brine L ow viscosity and 
freezing point permit use as a 
fow temperature cooling liquid. 


Most widely used low pressure 
propellant for noivfood use Does 
not provide adequate piessure alone, 
so IS almost entirely used m blends 
with "Freon" 12lc), 


Occasionally used as a solvent 
("Freon" MF) Blowing agent for 
foamed plastics’'^’ 


"F reon" 1 1 3 


fn commercial and industrial air 
conditioning and process water or 
brine chilling using centnlugal 
compressors, particularly in small 
tonnage applications. 


Solvent in some aerosol lormulations, 
usually propelled with "freon" 12 


E xlensively as a solvent ("Freon" TF) 
alone and »ri special purpose 
formulations for a wide range of 
critical cleaning needs In cutting 
fluid formulations. Valclene 
dry cleaning, etc. 



(a) A nurntM?f ot a/eoiropes l‘ Preon” 500, ' Prpon 507 etc ) are availdhte lor refrigeration use Bulletini describing ihe compostiion. 
propeMies and uses of these mixtures ere available on request. 

Ibl Normally aerosol propellants a»e hlemlerf to give the required vapor pressure end solubihry retiuiremen tt 

l c) "P reoo" It S IS a siatnli/ed grade trequently used m formulations containing alcohols Or waier 

ld) ‘‘Freoo 1 7 and ■' F reon ■ t1 are widely used as blowing agents for a range of toamed plasi ics m which ihey provide excellent cell structure 
and ihe trapped fluorocadron m closed celt loams signdicanily imoroves insulation and recovery properties ol the team 
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APPENDIX D 



THERMOCOUPLE CALIBRATION 



A. EQUIPMENT USED 

The equipment used in thermocouple calibration is shown 
in Figure D.l, and a brief description of each component is 
given below: 

1 . Thermocouple Wire 

Type-T ( copper- constantan) Teflon- coated wire of 
0.254 mm (0.01 in.) in diameter was used for all thermocou- 
ples . 



2 . Calibration Bath 

A Thermos flask was used as the calibration bath. In 
order to maintain an isothermal temperature distribution, a 
motor-driven mixer was used. Observations during the cali- 
bration procedure showed bath temperature fluctuations to be 
± 0.002 K. 

3 . Thermocouple Readout 

A Hewlett Packard 3497A automatic data acquisition/ 
control system and a Hewlett Packard 9826 computer were used 
to read, analyze and record the calibration data. 

4 . Reference Temperature 

A Hewlett Packard 2804A quartz thermometer was used 
to measure the bath temperature. This quartz thermometer had 
a resolution of 0.0001 K, while the manufacturer- guaranteed 
accuracy was better than ± 0.03 K. 
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B. PREPARATION FOR CALIBRATION 



1 . Thermocouple Preparation 

A total of five thermocouples were prepared for the 
calibration run. Two thermocouples were made from the 
beginning of a spool of copper- cons tantan wire; two were 
made from the end; and one was made at the mid section using 
the following procedure: 

• The Teflon insulation was removed for a length of about 
4 mm from one end of aim long piece of wire and a 
thermocouple bead was made using a Dynatech Corporation 
thermocouple welder. 

• The other end of the thermocouple wire was connected to 
the data acquisition and control unit through a junc- 
tion box. 

2 . Computer Program 

A short computer program (TCAL) was written to 
accept the thermocouple readings through the data acquisi- 
tion system and the bath temperature through the digital 
quartz thermometer. A listing of the computer program, TCAL 
is provided In Appendix G. This program prints all data and 
the discrepancy (i.e., the Quartz thermometer reading minus 
the thermocouple reading) as well as it stores the data on a 
computer disk. In order to convert the e.m.f. values to 
temperature, the manufacturer’s conversion equation (D.l) 
was used in the computer program TCAL: 



T = 2Lq+ + a3E' (D.l) 

+ a^E” + a^E^ + a^E®+ 
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where : 



T = temperature 
a<3= 0.100860910 
ai= 25727.94369 
a2 = -767345.8295 
a3= 78025595.81 
a4= -9247486589 
a5= 6.97688E+11 
a^= -2.66192E+13 
aj= 3.94078E+14 
E = thermocouple 



(°C) 



reading (volts) 



C. CALIBRATION PROCEDURE 

• Since the temperature measurements during this investi- 
gation ranged from 25 °C to 70 °C, thermocouple cali- 
bration was performed in this region. First, hot water 
(at about 75 °C) was added to the thermos flask cali- 

• bration bath. Following proper mixing, all thermocouple 
readings were recorded through the data acquisition 
system and the bath temperature was measured by the 
quartz thermometer. To obtain lower temperatures, 
small quantities of cold water were gradually added to 
the bath. 

• The discrepancy (i.e., the Quartz thermometer reading 
minus the thermocouple reading) was plotted against the 
thermocouple reading as shown in Figure D.2. The 
second-order polynomial curve (D.2) shown in this 
figure was generated using the data from all 5 
thermocouples . 
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DCP 



-6 , 7422934E-02 + 9 . 0277043E- 03 T 



(D.2) 



9 . 3259917E-05 T" 



where : 



T = thermocouple reading 

(from equation D.l) (°C) 
DCP = discrepancy (K) 



Thus, the corrected temperature values were 
using equation (D.3): 



obtained 



Tb = DCP + T 



(D.3) 



where : 



= actual calculated temperature (°C) 
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Figure D.l Schematic of the Calibration Devices. 
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Figure D.2 Thermocouple Calibration Curve. 



APPENDIX E 



DATA REDUCTION PROGRAM 



1000! f^ILE NttliE: DRPRl 
1010! DOTE; October 19. 198 A 

I0:0! REVISED: November 1. 1980 

UV?0! 

lOOO BEEP 

1050 PRINTER IS 1 

1060 PRINT USING "OX. ""Select option: 

1070 PRIfJT USING "GX.""0 Taking data or re-process i nq previous data""" 

lORO PRINT USING "bX.""1 Plotting data 

1090 INPUT I dp 

1100 IP ldp = 0 THEN COLL Mam 

1110 IF Idp-1 THEN CALL Plot 

1120 END 

1 1 30 SUB Main ■ 

1100 COM /Cc/ C(7) 

1150 DIM Emf ( 12) . T(I2) .DIa<2) .D2a<2) .Dia<2) .Doa(2) .La(2) .Lua(2) ,Kcua(2) 
1160 DATA 0. 10086091 .25727 .90369,-767305.8295.78025595.81 
1170 DATA -9207086589. 6. 97688E+1 1 .-2. 66192E+13. 3. 90078E+10 
1 180 READ C(“) 

I 190 PRINIER IS 701 
1200 CLEAR 709 
1210 BEEP 

1220 INPUT "ENTER MONTH. DATE AND TIME (MM ; DD :HH : MM :SS)" ,Date$ 

1230 OUTPUT 709:"TD" :DateS 
1200 OUTPUT 709:"TD" 

1250 ENTER 709:Date$ 

1260 PRINT 
1270 PRINT 
1280 PRIfIT 

1290 PRINT " Month, date and time :";Date'£ 

1300 PRINT 

nio PRINT USING "10X.""NOTE: Program name : DRP 

1320 BEEP 

1330 INPUT "ENTER DISK NUMBER”. Dn 

1300 PRINT USING "1GX,""Disk number =■ "".ZZ”:Dn 

1 350 PEEP 

1360 INPUT "ENTER INPUT MODE < 0 = 3050A . I =F ILE ) " . Im 
1370 IF lm=0 THEN 
1380 BEEP 

1390 INPLII "GIVE A NAME FOR THE RAN DATA FILE" . D2_f i leS 
1 000 PR INI USING "16X.""New file name: "" . 1 OA" : D2_f i leS 
1010 CREATE BOAT D2_file$,20 
1020 ASSIGN ®File2 TO D2_file$ 

H30! 

1000! DUMMY FILE UNTIL Nrun KNONN 
1060 D1_f i leS="DUMMY" 

1060 CREAIE BDAT D1_flle$,.?0 
1 070 ASSIGN fjiFllel TO D1_fiieS 
1080 OUTPUT ®F i lei :Date5 
1090 BEEP 

1500 INPUT "GIVE A NAME FOR THE PLOT FILE" .P_f i le$ 

1510 CREATE BDAT P_f i leS,5 
1520 ASSIGN iiPlot TO P_file$ 

1530 BEEP 

1500 INPUT "ENTER NUMBER OF DEFECTIVE ICS ( 0 = DEFAULT ) " . Id tc 
1550 IF Idtc’O THEN 
1560 Ldtcl=0 
1570 Ldtc2-0 

1580 PRINT USING "I6X,"”No defective TCs exist' 

1590 END IF 
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loOO IF Idtc=) THEN 
1610 BEEP 

1620 INPUT "ENTER DEFECTIVE TC LOCft I ION" . Ld tc 1 

1630 PRINT USING "16X,""TC is defective at location ”",D":Ldtcl 

1640 Ldtc2=0 

1650 END IF 

1660 IF Idtc=2 THEN 

1670 BEEP 

1680 INPUT "ENTER DEFECTIVE TC LOCAT IONS" . Id tc 1 . Ld tc2 

1680 PRINT USING "16X,""TC are defective at locations "".D.4X,D"; 

1700 END IF 

1 710 IF Idtc>2 THEM 

1720 BEEP 

1730 PRINTER IS 1 

1740 BEEP 

1750 PRINT "INVALID ENTRY" 

1760 PRINTER IS 701 
1770 GOTO 1530 
1780 END IP 

1790 OUTPUT ^Filel :LdtcI ,Ldtc2 
1800! Im= 1 option 
1810 ELSE 
1820 BEEP 

1830 INPUT "GIVE THE NAME OF THE EXISTING DATA FILE" ,D2_f i leS 
1840 PRINT USING "16X.""01d file name: "" . 1 4A" ; D2_f i leS 
1850 ASSIGN SFile2 TO D2_file$ 

1860 ENTER @File2:Nru.n 
1870 ENTER §)F i le2 ; Do ld$ 

1880 BEEP 

1890 INPUT "GIVE A NAME FOR PLOT F ILE" . P_f i leS 
1900 CREATE BDAT P_fileS.5 
1910 ASSIGN @Plot TO P_file'B 

1920 PRINT USING "I6X,""This data set taken on : . I 4A" : Do IdB 

1930 ENTER ®File2:Ldtc1 ,Ldtc2 

1940 PRINT USING " I 6X , ""Thermocoup 1 es were defective at locations 

Ld tc 1 . Ld tc2 

1950 ENTER ?iFile2:Itt 

1960 END IF 

1970 PRINTER IS I 

1980 IF lm=0 THEN 

1990 BEEP 

2000 PRINT USING "4X. ""Select tube type""" 

2010 PRINT USING "6X . ""0=Smooth 4 inch Ref 

2020 PRINT USING "GX . "" I =Smooth 4 inch soft Cu 

2030 PRINT USING "GX . ""2=Smoo th 8 inch soft Cu 

2040 INPUT Itt 

2050 IF Itt>2 THEM 
2060 BEEP 

2070 PRINT "INVALID ENTRY" 

2080 GOTO 2000 

2090 END IF 

2100 OUTPUT ®File1 :Itt 

2110 END IF 

2120 PRINTER IS 701 

2130 PRINT USING "1GX.""Tube Type is : "”.D":Itt 
2140 BEEP 

2150 INPUT "ENTER OUTPUT VERSION ( 0=LOMG . 1 -SHORT )". lov 
2160! 

2170! D1=Diameter at thermocouple positions 
2180 DATA .01 1 125. .01 1 1 125 . .01 1 43 
2190 READ Dla(‘) 



Ld tc 1 . Ld tc2 



:"".2(3D.4X)"i 
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2200 D1=D1a(Itt) 

2210 ! 

2220? D2=Diameter of test section to the base of fins 
2230 DATA . 01 5875 0 1 5875 .. 0 1 5875 
2240 READ D2a(-) 

2250 D2=D2a(Itt) 

2260? 

2270! Di=Inside diameter of unenhanced ends 
2280 DATA .0127. .0127, .0127 
2290 READ Dia(>*) 

2300 Di=Dia( Itt> 

2310! 

2320! Do=0utside diameter of unenhanced ends 
2330 DATA . 01 5875 .. 0 1 5875 .. 0 1 5875 
2340 READ Doa(-) 

2350 Do=Doa(Itt) 

2360! 

2370! L=Length of enhanced surface 
2380 DATA . 1016, . 1016. .2032 
2390 READ La(») 

2400 L=La(Itt) 

2410? 

2420! Lu=Length of unenhanced surface at the ends 
2430 DATA .0254, .0254. .0762 
2440 READ Lua(*> 

2450 Lu=Lua(Itt) 

2460! 

2470! Kcu'Thermal Conductivity of tube 
2480 DATA 398.344.344 
2490 READ Kcua(-) 

2500 Kcui=Kcua( 1 1 1 ) 

2510 A=PI-(Do 2-Di 2)/4 

2520 P=PI*Do 

2530 J=I 

2540 Sx=0 

2550 .Sy = 0 

2560 Sxs=0 

2570 Sxy=0 

2580 Repeat : ! 

2590 IF lm=0 THEN 

2600 ON KEY 0.15 RECOVER 2580 

2610 PRINTER IS 1 

2620 PRINT USING "4X. ""SELECT OPTION 

2630 PRINT USING "9X.""0=TAKE DATA 

2640 PRINT USING "9X,""1=SET HEAT FLUX 

2650 PRINT USING "9X.""2=SET Tsat 

2660 PRINT USING "4X,""N0TE: KEY 0 = ESCAPE 

2670 BEEP 

2680 INPUT Ido 

2690 IF ldo=0 THEN 3680 

2700! 

27!0! LOOP TO SET HEAT FLUX 

2720 IF Ido=1 THEN 

2730 OUTPUT 709;"AR AF62 AL63 VR5" 

2740 BEEP 

2750 INPUT "ENTER DESIRED 0dp",Dqdp 

2760 PRINT USING "4X , ""DESIRED Odp ACTUAL Odp 

2770 Err=1000 

2780 FOR 1=1 TO 2 

2790 OUTPUT 709: "AS SA" 

2800 Sum=0 



92 



2a 10 FOR Ji=I TO b 
2820 ENTER 709:E 
2830 Sum=Sum+E 
2840 NEXT Ji 

2850 IF I=I THEN Oolt=Sum*5 
2860 IF 1=2 THEN Amp=E 
2870 NEXT I 

2880 Aqdp=Volt»Amp/(PI‘D2«L) 

2890 IF ABS(Aqdp-Dqdp)>Err THEN 

2900 IF Aqdp>Dqdp THEN 

2910 BEEP 4000. .2 

2920 BEEP 4000.. 2 

2930 BEEP 4000. .2 

2940 ELSE 

2950 BEEP 250. .2 

2960 BEEP 250.. 2 

2970 BEEF 250. .2 

2980 END IF 

2990 PRINT USING ”4X . MZ . 3DE . 2X . HZ . 3DE" : Dqdp . Aqdp 

3000 WAIT 2 

3010 GOTO 2780 

3020 ELSE 

3030 BEEP 

3040 PRINT USING "4X . HZ . 3DE . 2X . HZ . 3DE" ; Dqdp . Aqdp 

3050 Err=500 

3060 WAIT 2 

3070 GOTO 2780 

3080 END IF 

3090 END IF 

3100! 

3! 10! LOOP TO SET Tsat 

3120 IF Ido=2 THEN 

3130 OUTPUT 709:"AR AF33 AL33 VR5" 

3140 OUTPUT 709:"AS SA" 

3150 BEEP 

3160 INPUT "ENTER DESIRED Tsaf.Dtld 

3170 PRINT USING "4X.’”'D Tsat A Tsat Diff Psat" 

3180 01d=0 

3190 ENTER 709:Eliq 

3200 Atld = FNTvsv<El iq) 

3210 Psat=FNPsat(Atld) 

3220 IF ABS(Atld-Dtld)>.2 THEN 

3230 IF Atld>Dtld THEN 

3240 BEEP 4000.. 2 

3250 BEEP 4000, .2 

3260 BEEP 4000.. 2 

3270 ELSE 

3280 BEEP 250. .2 

3290 BEEP 250.. 2 

3300 BEEP 250.. 2 

3310 END IF 

3320 Err=Atld-01d 

3330 01d=Atld 

3340 PRINT USING "4X . 4 < HDD . DD . 3X X" : D t Id . A t Id . Er r .Psat 
3350 WAIT 2 
3360 GOTO 3190 
3370 ELSE 

3380 IF ABS(Atld-Dtld)> . 1 THEN 
3390 IF Atld>Dtld THEM 
3400 BEEP 3000. .2 
3410 BEEP 3000. .2 
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3420 ELSE 
3430 BEEP 800, .2 
3440 BEEP 800.. 2 
3450 END IF 
3460 Err=8tld-01d 
3470 Old=Atld 

3480 PRINT USING "4X . 4( MDD . DD . 3X) " : D t Id . 8 t Id . Er r . Psat 

3490 WAIT 2 

3500 GOTO 3J90 

3510 ELSE 

3520 BEEP 

3530 Err=Atld-01d 

3540 01d=Atid 

3550 PRINT USING ”4X . 4 (MDD . DD . 3X ) " : Ot Id . At Id . Er r . Psat 

3560 WAIT 2 

3570 GOTO 3190 

3580 END IF 

3590 END IF 

3600 END IF 

3610! ERROR TRAP FOR Ido OUT OF BOUNDS 
3620 IF Ido>2 THEN 
3630 BEEP 
3640 GOTO 2620 
3650 END IF 
3660! 

3670! TAKE DATA IF lm=0 LOOP 
3680 BEEP 

3690 INPUT "ENTER BULK OIL X".Bop 
3700 OUTPUT 709:"AR AF25 AL36 VR5" 

3710 FOR I=I TO 12 
3720 OUTPUT 709: "AS SA" 

3730 Sum=0 
3740 FOR Ji=1 TO 20 
3750 ENTER 709:E 
3760 Sam=Sum+E 
3770 NEXT Ji 
3780 Emf (I)=Sum/20 
3790 NEXT I 

3300 OUTPUT 709:"AR AF62 AL63 VR5" 

3810 FOR 1=1 TO 2 
3820 OUTPUT 709: "AS SA" 

3830 Sum=0 
3840 FOR Ji=l TO 20 
3850 ENTER 709 :E 
3860 Sum=Sum+E 
3870 NEXT Ji 

3880 IF 1=1 THEN Ur=Sum/20 
3890 IF 1=2 THEM Ir=Sum/20 
3900 NEXT I 
3910 ELSE 

3920 ENTER ?F i le2 : Bop . To ld$ . Emf < • ) . Vr . I r 
3930 END IF 
3940 ! 

3950! CONVERT emf'S TO TEMP . VOLT . CURRENT 

3960 Tua=0 

3970 FOR I=I TO 12 

3980 IF Idtc>n THEN 

3990 IF I=Ldtcl OR I=Ldtc2 THEN 

4000 T(I)=-99.99 

4010 GOTO 4060 

4020 END IF 
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lNL) II- 

4040 T(I)=FNTvsv(Emf<I)) 

4050 IF I<9 THEN Twa=Tua+T(D 

4060 NEXT I 

4070 Tua=Tuia/(8* Idtc ) 

4080 Tld=T(9) 

4080 Tv=(T< 10)+T(1 1 ))/2 
4100 Tsump= T( I 2) 

4110 Amp= I r 
4120 Volt=Vr-25 
4130 0 = Vo 1 1 * Amp 

4140 IF ltt=0 THEN 
4150 Kcij.=FNKcu(Tua) 

4160 ELSE 

4170 Kcu.= Kcua< 1 1 1 > 

4180 END IF 
41 80! 

4200? FOURIER CONDUCIION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
4210 Tu. = Twa-0-L0G(D2/D1 ) / ( 2‘P I *Kca*L ) 

4220 The tab=Tu- T 1 d 
4230 IF Thetab<0 THEN 
4240 BEEP 

4250 INPUT •*TWALL<TSAT (0 = CONTINUE. i=END)'Mev 

4260 IF lev=0 THEN GOTO 2580 

4270 IF Iev=1 THEN 5030 

4280 END IF 

4280! 

4300! COHPUTE VARIOUS PROPERTIES 
4310 Tf ilm=FNTf ilm(T(j,Tld) 

4320 Rho=FNRho(Tfilm) 

4330 Mu=FNMu(Tf i Im) 

4340 K=FNK(Tfilm) 

4350 Cp=FNCp< Tf i Im) 

4360 Beta=FNBeta< Tf i Im) 

4370 Ni=Mu/Rho 
4380 AIpha=K/<Rho»Cp) 

4380 Pr=Ni/Alpha 
4400 Psat=FNPsat ( T 1 d ) 

4410! 

4420! COMPUTE HEAT TRANSFER COEFFICIENT 
4430 Hbar=180 

4440 Fe= < Hbar ‘P/ < Kca*A ) ) ' . 5*La 
4450 Tanh=FNTanh(Fe) 

4460 Theta=Thetab* Tanh/Fe 
4470! PRINTER IS 701 

4480! PRINT USING "4X , 7 < I X . MZ . 3DE ) " ; Hbar . Fe . Tanh . The tab . The ta . Be ta , N i 
4490 Xx=(9,81*Beta'‘Thetab*Do 3*Tanh/<Fe*Ni*Alpha)) .166667 
4500 Yy = ( I +< .559/Pr) '(9/16) ) (8/27) 

4510 Hbarc=K/Do‘( . 6+ . 387-X x / Yy ) 2 

4520 IF ABS((Hbar-Hbarc)/Hbarc)>.001 THEN 

4530 Hbar= ( Hbar +Hbarc ) • . 5 

4540 GOTO 4440 

4550 END IF 

4560 0 1 = ( Hbar*P»Kcu*‘A ) . 5-* The tab» Tanh 

4570 0c=0-2-01 

4580 As=PI‘D2«L 

4590 Odp=Oc/As 

4600 Htube=Qdp/The tab 

4610! 

4620! RECORD TIME OF DATA TAKING 
4630 IF lm=0 THEN 
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4fc40 OUTPUT 709; "TD" 

4650 ENTER 709;Told$ 

4660 END IF 
4670f 

^680! OUTPUT DATA TO PRINTER 
4690 PRINTER IS 701 
4700 PRINT 
4710 IF lov=0 THEN 

4720 PRINT USING "1 IX, ""Data Set Number = . DDD . 2X . ""Bulk Oil 7. = "".DD.D.5X.I 

4A" : J , Bop . Told? 

4730 PRINT 

4740 PRINT USING "1IX.""TC No : 1 2 3 4 5 6 7 

8 

4750 PRINT USING "1!X.""Tenp ; , 8 ( I X . HDD . DD ) " : T ( I ) . T < 2 ) . T ( 3 ) . T ( 4 ) . T ( 5 ) , T ( 6 ) . T ( 

7).T<8) 

4760 PRINT USING "!1X,"“ Twa Tliqd Thetab Tvapr Psat Tsumo 

4770 PRINT USING " 1 1 X . 2 ( MDD , DD . 1 X ) . MZ . DDDE . 1 X . 2 ( 1 X . HDD . DD ) . 2X . HDD . D" : Tua . T 1 d . Th 
e tab , Tv . Psat . T sump 

4780 PRINT USING "11X,"" Htube Qdp 

4790 PRINT USING "1 1 X , MZ . 3DE . 1 X . MZ . 3DE" : Htube . Qdp 

4800 END IF 

4810 IF Iov=1 THEN 

4820 IF J=l THEN 

4830 PRINT USING "1 IX, ""RUN No OilX Tsat Htube Odp Thetab 

4840 END IF 

4850 PRINT USING " 1 2X . DDD . 4X , DD , 2X . NDD . DD . 3d X , MZ . DDDE> " : J , Bop , Tld . Ht LLbe , Odp . Th 
e tab 

4860 END IF 
4870 IF lm=0 THEM 
4880 BEEP 

4890 INPUT "OK TO STORE THIS DATA SET ( ! =Y ,0=N) ?" ,0k 
4900 END IF 

4910 IF 0k=l OR Im=l THEN J=J+1 

4920 IF Okd AND lm=0 THEM OUTPUT ®F i le 1 : Bop . To ld$ , Emf ( ‘ ) . Ur , Ir 
4930 IF Im=1 OR Ok = 1 THEM OUTPUT ®P lot . -Qdp . Thetab 
4940 IF lm=0 THEM 
4950 BEEP 

4960 INPUT "WILL THERE BE ANOTHER RUN d = Y . 0 = N ) . Go_on 
4970 Nrun=J 

4980 IF Go_on01 THEN 5030 
4990 IF Go_on= 1 THEN Repeat 
5000 ELSE 

5010 IF J<Nrun+1 THEM Repeat 

5020 END IF 

5030 IF lm=0 THEN 

5040 BEEP 

5050 PRINT 

5060 PRINT USING " 1 0X , ""NOTE : "".ZZ."" data runs uere stored in file "".I0A":J- 
1 ,D2_fileT. 

5070 ASSIGN 9Filel TO * 

5080 OUTPUT ®Fi le2;Nrun-1 

5090 ASSIGN ?FileI TO D1 file? 

5100 ENTER ®F i lei : Date? . Ldtc 1 , Ldtc2 , It t 
5110 OUTPUT 9File2:Date?,Ldtcl ,Ldtc2. Itt 
5120 FOR 1=1 TO Nrun-1 

5130 ENTER 9F i le I : Bop , To Id? , Emf < • ) , Ur , I r 

5140 OUTPUT ®File2:Bop.Told?,Emf (•) ,Ur,Ir 

5150 NEXT I 

5160 ASSIGN AFilel TO 

5170 PURGE "DUMMY" 

5180 END IF 
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-> I JU DLC.f 

:d^^UU PRlNi 

5210 PRINT USING "10X.""NOTE: X-Y pairs uere stored in plot data file 

.P_f ileT 

5220 ASSIGN ?File2 TO * 

5230 ASSIGN @Plot TO • 

5240 BEEP 

5250 INPUT "LIKE TO PLOT DATA ( 1 =Y , 0=N)?" ,0k 

5260 IF 0k=1 THEN 

5270 CALL Plot 

5280 END IF 

5290 SUBEND 

5300! 

5310! CURVE FITS OF PROPERTY FUNCTIONS 
5320 DEF FNKcu.(T) 

5330! OFHC COPPER 250 TO 300 K 

5340 Tk=T+273.15 !C TO K 

5350 K=434- . I l2»Tk 

5360 RETURN K 

5370 FMEND 

5380 DEF FNMu.(T> 

5390! 170 TO 360 K CURVE FIT OF VISCOUSITY 

5400 Tk=T+273.15 !C TO K 

5410 Hu=EXP(-4.4636+( !011.47/Tk))^l .OE-3 

5420 RETURN Mu 

5430 FNEND 

5440 DEF FNCp<T) 

5450! 180 TO 400 K CURVE FIT OF Cp 

5460 Tk=T+273.15 !C TO K 

5470 Cp= .40188+1 .65007E-3*Tk +1 .51494E-6-Tk 2-6 . 67853E- 1 0«Tk '3 

5480 RETURN Cp 

5490 FNEND 

5500 DEF FNRho(T) 

5510 Tk=T+273. 15 !C TO K 

5520 X=l-(1 .8»Tk/753.95) !K TO R 

5530 Ro=36. 32 + 61 . I 464 1 4*X ' < 1 73 ) + 1 6 . 4 1 80 1 5-X+ I 7 . 476838*X .5+1.1 I9828*X'? 

5540 RETURN Ro 
5550 FNEND 
5560 DEF FMPr(T) 

5570 Pr = FNCpa)»FNMu(T)/FNK(T) 

5580 RETURN Pr 
5590 FNEND 
5600 DEF FNK(T) 

5510! T<360 K WITH T IN C 
5620 K= .071 - .000261 *T 
5630 RETURN K 
5640 FNEND 
5650 DEF FNTanh(X) 

5660 P=EXP(X> 

5670 Q=EXP(-X) 

5680 Tanh=(P-0)/(P+Q) 

5690 RETURN Tanh 

5700 FNEND 

5710 DEF FNTvsv(V) 

5720 COM /Cc/ C(7) 

5730 T=C(0) 

5740 FOR 1=1 TO 7 
5750 T=T+C<I)*V I 
5760 NEXT I 

5 770 T = T-6. 7422934E-2 + T‘(9.0277043E-3-TM-9.3259917E-5)) 
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5780 RETURN T 

5790 FNEMD 

5800 DEF FMBetaCT) 

5810 Rop=FMRho(T+. I) 

5820 Rom=FNRho< T- . 1 ) 

5830 Beta=-2/(Rop+Rom)»(Rop-Rom)/.2 
5840 RETURN Beta 
5850 FNEMD 

5880 DEF FNTf 1 lm(Tu. Tld) 

5870 Tf i lm=<Tu+Tld)/2 
5880 RETURN Tfilm 
5890 FNEMD 
5900 DEF FNPsat(Tc) 

5910? 0 TO 80 deg F CURVE FIT OF Psat, 

5920 Tf=1.8*Tc+32 

5930 Pa=5.945525 + Tf *( . 1 5352082 + Tf « ( I . 4840963E-3+ Tf **9 . 6 1 5067 1 E-6 ) ) 
5940 Pg=Pa-14.7 

5950 IF Pq>0 THEN ! +=PSIG.-=in Hg 

5960 Psat=Pg 
5970 ELSE 

5980 Psat=Pq*29.92/l4.7 

5990 END IF’ 

6000 RETURN Psat 
6010 FNEND 
6020 SUB Plot 
6030 DIM C(9> 

6040 INTEGER li 
6050 PRINTER IS 1 



6060 

6070 

6080 


Idv= 1 

BEEP 

INPUT 


"LIKE DEFAULT 


VALUES FOR PLOT ( l = Y.0 = 


6090 

6100 


BEEP 

PRINT 


USING "4X." 


"Se 


lect Uption: 


6110 


PRINT 


USING "4X."' 


"0 


q versus delta-T 


6120 


PRINT 


USING "4X." 


"1 


h versus deita-T 


6130 


PRINT 


USING "4X." 


"2 


h versus q 


6140 


PRINT 


USING "4X," 


"3 


h-ratio versus delta-T 


6150 


INPUT 


Opo 







6160 IF 0po=3 THEN 
6170 BEEP 

6180 INPUT "SELECT TUBE DIAMETER (0=. 75. 1 = 1.0 IN)'Mtd 

6190 END IF 

6200 PRINTER IS 705 

6210 IF IdvOI THEM 

6220 BEEP 

6230 INPUT "ENTER NUMBER OF CYCLES FOR X-AXIS". Cx 
6240 BEEP 

6250 INPUT "ENTER NUMBER OF CYCLES FOR Y-AXIS". Cy 
6260 BEEP 

6270 INPUT "ENTER MIN X-VALUE (MULTIPLE OF 10)",Xnun 
6280 BEEP 

6290 INPUT "ENTER MIN Y-VALUE (MULTIPLE OF I0)".Yt:iin 

6300 ELSE 

6 510 Cy = 3 

6320 IF 0po=0 them 

6330 Cx=2 

6340 Xm 1 n= 1 

6350 Ynun=1000 

6360 END IF 

6370 IF 0po=1 THEN 

6380 Cx=2 
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DJdU 

6400 

6410 

6420 

6430 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

66 1 0 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

6810 

6820 

6830 

6840 

6850 

6860 

6870 

6880 

6890 

6900 

6910 

6920 

6930 

6940 

6950 

6960 

6970 

6930 



Xm 1 n= I 
Ymin= 1 00 
END IF 

IF 0po=2 THEN 
Cx = 3 

Xnu n= 1 000 
Ymin= 1 00 
END IF 

IF 0po=3 THEN 

Cx = 2 

Cy=3 

Xm i n= 1 

Y m 1 n = 1 

END IF 

END IF 

BEEP 

PRINT "IN:SP1:IP 2300,1800.8300.6800;" 

PRINT "SC 0,100.0.100:TL 2.0:" 

Sfx=l00/Cx 

Sfy=100/Cy 

PRINT ”PU 0.0 PD” 

Nn = 9 

FOR 1=1 TO Cx+1 
Xat=Xmin*1 0' ( I- 1 > 

IF I=Cx+l THEN Nn=1 
FOR J=1 TO Nn 

IF J=1 THEN PRINT "TL 2 0" 

IF J=2 THEN PRINT "TL 1 0" 

Xa=Xat«J 

X=LGT(Xa/Xmin)»Sf X 
PRINT "PA":X.",0: XT:" 

NEXT J • 

NEXT I 

PRINT "PA I00.0;PU:" 

PRINT "PU PA 0.0 PD" 

Nn = 9 

FOR 1=1 TO Cy+1 
Yat=Ymin»10' (I-l ) 

IF I=Cy+l THEN Nn= 1 
FOR J=1 TO Nn 

IF J=l THEN PRINT "TL 2 0" 

IF J=2 THEN PRINT "TL 1 0" 

Ya= Yat • J 

Y=LGT( Ya/Ymin)-Sfy 
PRINT "PA 0,":Y."YT" 

NEXT J 
NEXT I 

PRINT "PA 0.100 TL 0 2" 

Nn = 9 

FOR 1=1 TO Cx+1 
Xat = Xm 1 n* 1 0 ( I - 1 ) 

IF I=Cx+l THEN Nn=1 
FOR J=1 TO Nn 

IF J=1 THEN PRINT "TL 0 2" 

IF J>1 THEN PRINT "TL 0 1" 

Xa=Xat • J 

X=LGT(Xa/Xmin)*Sf X 
PRINT "PA":X.". 100; XT" 

NEXT J 
NEXT I 
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6990 PRINT "Pfl 100.100 PU PA !00.0 PD” 

7000 Nn=9 

7010 FOR 1=1 TO Cy+1 
7020 Yat=Ymin>10' < 1-1) 

7030 IF I = Cy + 1 THEN Nn= I 

7040 FOR J=1 TO Nn 

7050 IF J=I THEN PRINT "TL 0 2" 

7060 IF J>1 THEN PRINT "TL 0 1" 

7070 Ya=Yat*J 

7080 Y=LGT(Ya/Ymin)*Sfy 

7090 PRINT "PD PA 100.".Y."YT" 

7100 NEXT J 

7110 NEXT I 

7120 PRINT "PA 100. 100 PU" 

7130 PRINT "PA 0.-2 SR 1 .5.2" 

7140 Ii=LGT(Xmin) 

7150 FOR 1=1 TO Cx+1 
7160 Xa=Xmin* I 0 ' ( I- 1 ) 

7170 X=LGT(Xa/Xmin)*Sf X 
7180 PRINT "PA”:X.".0;" 

7190 IF li> = 0 THEN PRINT "CP -2 . -2 : LB 1 0 : PR -2 . 2 : LB" : I i : 

7200 IF IKO THEN PRINT "CP -2 . -2 :LB1 0 :PR 0 . 2 : LB" : 1 1 : 

7210 Ii=Ii+1 

7220 NEXT I 

7230 PRINT "PU PA 0.0" 

7240 Ii=LGT(Ymin) 

7250 Y10=10 

7260 FOR 1=1 TO Cy+1 

7270 Ya=Y(iun*10 (1-1) 

7280 Y=LGT(Ya/Ymin)*Sfy 
7290 PRINT "PA 0.":Y."" 

73JOO PRINT "CP -4.-.25:LB10:PR -2 . 2 : LB" : 1 1 : 

7310 Ii=Ii+l 
7320 NEXT I 
7330 IF IdvOI THEN 
7340 BEEP 

7350 INPUT "ENTER X-LABEL" . X labels 
7360 BEEP 

7370 INPUT "ENTER Y-LABEL" . Y labels 

7380 ELSE 

7390 IF 0po=0 THEN 

7400 Xlabel$=”Tu-Tsat (K)" 

7410 Ylabel$="q (W/m'2)" 

7420 END IF 

7430 IF Opo=l THEN 

7440 Xlabel$="Tu-Tsat <K)" 

7450 Ylabel$="h (W/m'2. K)" 

7460 END IF 

7470 IF 0po=2 THEN 

7480 Xlabel$="q (W/m'2)" 

7490 YlabelS="h (W/m 2.K)" 

7500 END IF 

7510 IF l.lpo=3 THEN 

7520 XlabelS = "Tii.-Tsat (K)" 

7530 YlabelS=”h(enh)/h(smooth)" 

7540 END IF 
7550 END IF 

7560 PRINT "SR 1.5.2:PU PA 50.-16 CP” ; -LEN( X labe IS ) /2 : "0 : LB" : X labe IS 
7570 PRINT "PA -14.50 CP 0 . " : -LEN( Y labe IS )/2*5/6 : "DI 0 . 1 ;LB" ; Y labe IS 
7580 PRINT "CP 0.0 DI" 

7590 Repeat : ? 
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/bUU 

7610 

7620 

7630 

7540 

7650 

7660 

7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7 770 

7780 

7730 

7800 

7810 

7820 

7830 

7840 

7850 

7860 

7870 

7880 

7890 

7900 

7910 

7920 

7930 

7940 

7950 

7960 

7970 

7980 

7990 

8000 

8010 

3020 

8030 

3040 

8050 

3060 

8070 

8080 

8090 

8100 

8110 

8120 

3130 

8140 

3150 

8160 

8170 

8180 

8190 



Bttr 

INPUT "UlftMT TO 
IF 0I< = 1 THEM 
BEEP 

INPUT "ENTER 
8SSIGN @File 
BEEP 
BEEP 

INPUT "ENTER 
BEEP 

INPUT "ENTER 
BEEP 

PRINTER IS I 
PRINT USING 
USING 
USING 
USING 
USING 
USING 



PLOT D8T8 FROM 8 FILE ( I = Y . 0=N) 7" . Ok 



THE NAME OF 
TO D files 



THE D8T8 FILE".D fileS 



THE 

THE 



BEGINNING 
NUMBER OF 



RUN NUMBER". Md 

X-Y P8IRS STORED" .Npairs 



"4X , 

”4X 

”4X, 

"4X 

"4X, 

"4X, 



"Se 

"1 

"3 

"5 

"6 

"7 



Sym 
IS 



PRINT 
=RINT 
^RINT 
^RINT 
^RINT 
INPUT 

PRINTER IS 705 
"RINT "PU DI" 

IF Sym=1 THEN PRINT "SM»" 

IF Sym=2 THEN PRINT "SM+" 

IF Sym=3 THEN PRINT "SMo" 

[F Md>I THEN 
•OR 1=1 TO (Md-1) 

INTER ?)File:Ya.Xa 
TEXT I 

:nd if 

•OR 1=1 TO Npairs 
INTER s)File:Ya.Xa 
IF Opo=1 THEN Ya=Ya/Xa 
[F 0po=2 THEN 
5 = Ya 

('a= Ya/Xa 
<a=0 



ec t a symbo 1 : 

Star 2 Plus sign'" 

Circle 4 Square 

Rombus 

R ight-side-up triangle' 
Up-side-doun triangle"' 



END IF 

IF 0po=3 THEN Ya= Ya/FNHsmooth( Xa. I td) 
X=LGT(Xa/Xmin)*Sfx 



Y=LGT( Ya/Ymin)*Sfy 
IF Sytn>3 THEN PRINT 
IF Sym<4 THEN PRINT 
PRINT "P8".X.Y."" 

IF Sym>3 THEN PRINT 
IF Sym=4 THEN PRINT 
IF Sym=5 THEN PRINT 
IF Sym=6 THEN PRINT 
IF Sym=7 THEN PRINT 
NEXT 1 
BEEP 



"SM" 

"SR 1 .4.2.4" 

"SR 1 7 I S" 

"UC2,4'^99! 0.-8. -4. 0.0. 8.4.0:" 
"UC3. 0.99. -3. -6. -3. 6. 3. 6. 3. -6:" 
"UCO. 5. 3. 99. 3. -8. -6. 0.3.8:" 
"UCO, -5. 3. 99. -3. 8. 6. 0.-3. -8:" 



ASSIGN ?)File TO * 



GOTO 7600 
END IF 

PRINT "PU SM" 

BEEP 

INPUT "WANT TO PLOT A POLYNOMIAL ( I = Y . 0 = N ) 7" . Go_.:m 

IF Go_on=1 THEN 

BEEP 

INPUT "ENTER LOWER AND UPPER X-LIMI TS" . XI 1 . X 1 u 
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8200 

3210 

8220 

8230 

8240 

8250 

8260 

8270 

8280 

8290 

8300 

3310 

8320 

8330 

3340 

8350 

8360 

8370 

8380 

8390 

8400 

3410 

8420 

8430 

8440 

8450 

8460 

8470 

8480 

8490 

8500 

8510 

8520 

8530 

3540 

8550 

3560 



FOR Xx=0 TO Cx STEP Cx/200 
Xa=Xm 1 n"* 1 0 ' Xx 

IF Xa<Xll OR Xa>Xlu THEN 8290 
Ya=FNPoly(Xa) 

Y=LGT( Ya/Ymin)*Sfy 
X=LGT<Xa/Xmin)»Sfx 
IF Y<0 THEN Y=0 
IF Y>100 THEN GOTO 8290 
PRINT •*PPI'*,X. Y."PD" 

NEXT Xx 
END IF 

PRINT "PIJ PH 0.0 SPO” 

SUBEND 

DEF FNHsmoothCX . ltd) 

Hs=FNPoly (X)/X 

IF Itd=1 THEN Hs=Hs*. 83347 

RETURN Hs 

FNEMD 

DEF FMPoly(X) 

Poly=-4.4123718E+2-X*(6 .8123917E+2-X*3.7416863E+2) 

RETURN Poly 
FNEND 

DEF FNPvst( Tsteam) 

DIM K(8) 

DOTH -7.691234564.-26.08023696. -168. 1706546.64.23285504,-1 18.9646225 
DATA 4.16711732.20.9750676.1 .E9.6 
READ K(^) 

T=(Tsteam+273. 15)7647.3 
Sum = 0 

FOR N=0 TO 4 
Sum=Sum+K(N)*‘( 1-T)-'(N+1 ) 

NEXT N 

Br = Sum/( T'‘< 1+K(5)*‘(1-T)+K(6)*(I-T) 2))-(1-T)/(K(7)»(1-T) 2+K(8)) 
Pr=EXP(Br) 

P=22120000*Pr 
RETURN P 
FNEND 
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APPENDIX F 



AN EXAMPLE OF REPRESENTATIVE DATA RUN 



Month, date and time : I 2 : 03 : I 6 ; Ml ; 2^^ 



NOTE: Program name : DRP 
PisT number - 02 

Old file name: NH05 

Thi? data set taken on ; ll:02:IO:3fj:2P 
Thermocouples were defective at locations: 0 D 

T ube Type is : 2 



['ala Set Number = I Bulk Oil 7. ’ 0.0 

TC No; 12 ^ /I s 

Temp : 70.05 55.07 bl.04 70. ?3 7f!.P8 

Ima find Thetab Tvapr Peat 

Sb.30 10.03 5.560E+0I _I0.0I 3.86 

Htube Qdp 

i.505Er03 8.588E^0^ 

Oat a Set Number • 2 Bulk Oil 7. = 0.0 

TC No: 1 2 3 4 5 

Temp : 52.64 42.70 45.49 ^9- 95 58.26 

Tii/a Tliod Thetab Tvapr Psat 

49.84 10.07 3.939E>01 10.30 3.88 

Htube Qdp 

.I.I52EH13 4.538E + 04 

Data Set Number = 3 Bulk Oil 7- 0.0 

TC No: 1 2 3 ^ 5 

Temp : 34.69 30.25 31.40 33.49 37.59 

Tiiia Tliqd Thetab I\'apr Psat 

33.72 10.00 2.357E+01 10.69 3.83 

Htube Qdp 

7.420EH12 1.749E + 04 

Data Set Number = 4 Bulk 0 i 1 7. = 0.0 

TC Mo: 12345 

Temp : 27.30 24.92 25.40 26.60 29.04 

Itiia Tliqd Thetab Tvapr Psat 

26.78 10.07 1.663E^0I 11.17 3.88 

Htube Qdp 

5.2ME + 02 .8.6B4E + 03 

Data Set Number * 5 Bulk Oil 7. * 0.0 

TC No: I 2 3 'i 5 

Temp : 22 . 1 0 21.03 21.06 21.41 23.11 

Ttiia Tliqd Thetab Tvapr Psat 

21.88 9.95 1.189E+01 11.47 3.80 

Htube Qdp 

3.094E+02 3.680E+03 

Data Set Number • 6 Bulk 0 1 1 7. ' 0.0 



11:02:11:01:46 

6 7 P 

60.22 76.21 59.24 
1 sunip 

-13.0 



11:02:11 :30 ;52 

6 7 f: 

44.85 57.78 47,04 

Tsump 
- 12.0 



11:02:11: 38:07 

6 7 8 

31.38 38.00 32.96 

Tsump 

-11.9 



11:02:11 :49:56 

6 7 8 

25.25 29.30 26.40 

Tsump 

- II .9 



I ' :02: 1 I :57:02 

6 7 8 

21.55 23.11 21,65 

Tsump 
- 12 . 1 



11:02:12:05:46 
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TC No: 1 2 ^ S P, 7 P. 

Temp ; 18. M 17,88 17.82 17.87 19.16 18.03 18.03 16.80 
Twa Tliqd Thetab Tvapr Psat Tsamp 

17.95 10.19 7.737E+00 12.00 3.96 -12.2 

Htube Odp 

I.728E+02 I.337E+03 

Data Set Number = 7 Bulk 0 i 1 7. = 0.0 1 ! ; 02 : ! 2 : ! 0 : ! 8 

TC No: 1 2 3 0 8 6 7 8 

Temp : 18,58 18.08 18.10 10.87 18.66 18.00 18.06 10.58 

Tua Tliqd Thetab Tvapr Psat Tsump 

18.17 10.12 5.0O2E + 00 12.30 ,3.92 -12.2 

Htube Odp 

1.221E+02 6.186E+02 

NOTE; 07 X-Y pairs were stored in plot data file P08 
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APPENDIX G 



LISTING OF CALIBRATION COMPUTER PROGPJUI (TCAL) 



100 ? FILE NttME: TCAL 

no ! REVISED; Octobor d, 198^ 

1 LO ! 

130 COM /Cc/ C(7) 

100 DIM EmF(O) ,T(0) ,D(0) 

ISO DAIA 0. 10086091 . 25 73 7 . 90369 . - 76 7305 . 8295 . 7803SSA5 . 8 1 
160 DATA -92^7086589. 6. 97688E1 1 , -2. 66192E13. 3. 90078EM 
I 70 READ C(-) 

130 CLEAR 709 

190 BEEP 

200 IMPl.lT "ENTER MUNTH. DATE AND TIME < MM : DD : HH : MM : SS) " . P$ 
210 J = p 

220 OUTPUT 709:"TD":BS 

230 OUTPUT 709:"TD" 

2^0 ENTER 709:A$ 

24 1 BEEP 

242 INPUT "WANT A HARD COPY ( 1 = Y . O'TDP’Mhp 

244 IF Ihp=1 THEN PRINTER IS 701 

245 F lbp=0 THEN PRINTER IS I 

250 PRINT USING "1 OX , ""Month . date and time = "".MA-iAS 
260 BEEP 

270 INPUT "ENTER INPUT MODE (I-3054A. 2=FILE)".Im 

280 IF Im=1 THEN 

290 BEEP 

300 INPUT "GIVE A NAME FOR DATA F ILE" . D_f i le$ 

310 CREATE BDAT D_fUe$,20 
320 ELSE 

330 BEEP 

340 INPUT "GIVE NAME OF EXISTING F ILE" . D_f i I eS 
350 BEEP 

360 INPUT "ENTER NUMBER OF DATA RUNS STORED". Mrun 

361 BEEP 

362 INPUT "GIVE A NAME FOR OUTPUT FILE".Ofile$ 

363 CREATE BDAT Of i 1 p$. 10 

364 ASSIGN 9Fi leo TO Of l le$ 

365 K=1 

370 END IF 

380 ASSIGN j le TO D_fi le$ 

390 IF Iin=1 THEN 

400 BEEP 

410 ! INPUT "ENTER BATH TEMPERATURE" . T_bath 

411 OUTPUT 713:"T2R2E" 

412 WAIT 2 

413 ENTER 7l3:T_batb 

420 OUTPUT 709:"AR AF29 AL33 VRI" 

430 FOR 1=0 TO 4 

440 OUTPUT 709;"AS SA" 

450 ENTER 709:Emr<I) 

460 NEXT I 

470 ELSE 

400 ENTER i'F i 1 e : T_bat h . Emf ( ^ ) 

490 END IF 

500 j = jn 

510 FOR 1=0 TO 4 

520 T ( I )=FNTvsv<ABS(Emf ( I ) ) ) 

530 D(I)=T_bath-T( I) 

540 NEXT I 

541 IF K=1 THEM 

542 EOR L=D TO 4 

544 OUTPUT «Fi leo;T(L) .D(D 
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5^15 NEXT L 

548 GOTO 700 

549 END IF 

550 PRINT 

5G0 PRINf USING ” I OX , ""tia th temperature ‘ . 3D . 3D . (Deq O ;T_bath 

570 PRINT USING " 1 OX . The rmocoup le readings < Deg O: 

580 PRINT USING " I OX , 5 < 3D . DD , 3X ) . 35X” : T < • ) 

590 PRINT USING " 1 OX , " "D i sc repanc i es <Deg C): 

600 PRINT USING ” 11 X . 5 ( NZ . DD , 4X ) . 24X" : D ( » ) 

610 IF Im=2 THEN 700 
620 BEEP 

630 INPUT -riK TO UCCEPT THIS SET ( I =Y . 0 = N)'7‘* .Qks 

640 IF 0Fs=0 THEN 

650 J=J-1 

660 GUTn 390 

670 ELSE 

680 OUTPUT s>File:T_bath.Emf(*) 

690 END IF 

700 BEEP 

710 IF Im=l THEM 

720 INPUT "WILL THERE BE HHOTHER SET < 1 = Y . 0 =N ) V" . Go_on 
730 IF Go_on=1 THEM 390 

740 ELSE 

750 IF J<Nrun THEN 390 
760 END IF 

770 PRINT 

780 IF Im=1 THEN 

790 PRINT USING "I0X,""NOTE: "".DD."" data sets are stored in ^ i 1 e "".I4A";J,D 
_T i Ie$ 

800 ELSE 

810 PRINT USlNG " 1 0X , ""NOTE : Above analysis uas perTormed from file "".14A":D_ 
f i le$ 

811 PRINT USING "lOX,"" Output data are stored in Tile . I 4A" ; Of i 1 eS 

820 END IF 

•330 ASSIGN SFile TO • 

840 END 

850 DEF FNTvsv(Emf) 

860 COM /Cc/ C(7) 

870 Sum-C<0) 

8B0 FOR 1-1 TO 7 

390 Sum = Sum iC ( I ) '•Em f I 

900 NEXT I 

'310 RETURN Sum 

920 FNEND 
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APPENDIX H 



SAMPLE CALCULATION 

Data run number 5 (saturation temperature was 10 °C and 
heat flux was about 20 kW/m^) was chosen for the sample 
calculation . 

A. TEST-SECTION DIMENSIONS 



^^o 


= 0.01588 


(m) 


Di 


= 0.01270 


(m) 


Di 


= 0.01143 


(m) 


^2 


= 0.01588 


(m) 


L 


= 0.20320 


(m) 


Lu 


= 0.07620 


(m) 



B . MEASURED PARAMETERS 



Vs 


— 


3 . 58 


(volts ) 


Is 




2.27 


(volts ) 


Tl 


= 


34.69 


(°c) 


l2 




30.25 


(°c) 


l3 




31.40 


(°c) 


l4 




33.49 


(°c) 


Ts 




37.59 


(°c) 


l6 




31.38 


(°c) 


T7 


- 


38.00 


(°c) 


Is 


= 


32.96 


(°C) 


I s at 


r 


10.00 


(°C) 


kc 


= 


344.00 


(W/m.K) 



107 



C. OUTER WALL TEMPERATURE OF THE BOILING TUBE 



p = = 3.1416 X 0.01588 = 0.0499 (m) 

^ (D"o - D\)/4 = 3.1416 x (0.01588" - 0.0127")/4 

A ^ = 7 14 X 10 "" (m" ) 

V = 25 V 3 = 25 X 3.58 = 89.50 (volts) 

I = 15 = 2.27 ( amp ) 

Q = VI = 89.50 X 2.27 = 203.17 (W) 



T = T - Q,.( ln(D^/ Dj/2 ttr k ) 
wo avg ^ 2 1' c ' 

33.72 - 203. 17(ln(15. 88/11. 43)/2x3. 1416x203. 2x344) 
%o= 33.57 (°C) 

“ ^wo 

0 = Tb- T = 3.3.57 - 10.00 = 23.57 (°C) 

^ sat 

D. PROPERTIES OF R-114 AT FILM TEMPERATURE {REF. 24}, (REF. 



T. = (f + T J/2 = 21.78 (°C) = 294.78 (K) 

^ WO S 3. t 

log y (10 -"Ns/m" ) = -4.4636 + 1011.47/Tf 
y = exp(-4.4636 + (1011.47/294.78)) x 10 
y = 356 X 10 (Ns/m" ) 

P = 36.32 + 61.14 j"'" + 16.42 j + 17.48 j"'" + 

1.12 j" (Ib/ft") 



8 




25} 



j = 1 - /T 

T^ = Critical Temperature = 753,95 (R) 
T^ = 294.78 (K) = 530.6 (R) 

j = 1 - 530.60/753.95 = 0.296 
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P = 91.55 (Ib/ft^) = 1466.49 (kg/m^) 

V = u / P = 2.428 X 10 (mVs) 

k (W/m.K) = 0.0710 - 0.000261 Tf ("% in °C) 

k = 0.0653 (W/m.K) 

Cp(kJ/kg.K) = 0.4 + 1.65 x 10-'Tf+ 1.51 x 10 - 

6.68 X 10 (T^ in K) 

= 1002 (J/kg.K) 

0^ = k/ P = 4.444 X 10 "® (m^/s) 

3 = -(Ap/AT)/P 

3 = -(1465.95 - 1466.49) / ((294.98 - 294.78) x 1466) 

3 = 0.00184 (1/K) 

Pr = V / a = 5.464 

E. HEAT- FLUX CALCULATION 



Average natural- convection heat- transfer coefficient 
at non-boiling ends; 



h 






^ / 0 . 60 + 0.387 
o 



_ 1/2 

gcD^O Tjnht ./ v ■) Lu^*^ 

^ o B { K 



c'^c 



I 



, ( I 1/2 



c c 



[1 + (0 . 559/Pr) 






/ 



h = 282 (W/in'.K) 

Heat-transfer rate through non-boiling ends: 
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Q = (h p 0j^Tanh(mL^) 

m = (h p / k^A^)^ = 23.94 (1/m) 

= 13.15 (W) 

Q^os3 ^ ^ ^ ^ 13.15 = 26.30 (W) 

Heat flux through active boiling surface: 

Q = Q^- 203.17 - 26.30 = 176.87 (W) 

q = Q / A j3 

A, = D^L = 10.137 X 10-' (m") 
q = 176.87 / 0.010137 = 17447.96 (W/m^) 
h = q / 0j^ = 740.26 (W/m".K) 

The following are the results obtained from the computer 
by running the data reduction program (See Appendix F). 

q = 17490 (W/m^) 

0j^ = 23.57 (°C) 

h = 742 (W/m".K) 
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APPENDIX I 



UNCERTAINTY ANALYSIS 

The same data set (run number 5) that was used for the 
sample calculation was chosen for the uncertainty analysis; 
therefore, the dimensions of the test section, and the meas- 
ured and calculated parameters found in the sample calcula- 
tion were used in this analysis. All uncertainties are 
presented as a percentage of the calculated parameter. 

A. UNCERTAINTY IN SOURCE HEAT-TRANSFER RATE 
Qh = VI (W) 

I = I^ = 2.27 (amp) *51 = ±0.025 (amp) 

V = 25V ; = 3.58 (volts) 6Vo = ±0.05 (volts) 

s s ^ 

where : 

6 = uncertainty in measurement and calculation 

«Qh/Qh = 

SQh/Qh = ((0.05/3.58)^ + (0. 025/2. 27)U‘ 

= 1.80 percent 

B. UNCERTAINTY IN SURFACE AREA 
Ab = ttD^L 

Dq = 15.88 (mm) 5D^ = 0.1 (mm) 

L = 203.20 ( mm ) 6L = 0 . 1 ( mm ) 

6Ab/Ab = ((5DJD^)^ ^ (6L/L)^)^'^ 

6Ab/Ab = ((0.1/15.88)=" + (0.1/203.2)")"'=" 

= 0.63 percent 
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C. UNCERTAINTY IN WALL SUPERHEAT 



AT = 

T = 
sat 

XfO “ 

T = 
a vg 



5T 



X; o " '^s a t 

10.00 (°C) 
tvg- % (ln(D2 /Di ) / 2, L k^) 

( / « 



0.5 (°C) 



n=l 



where : 

= thermocouple readings (See Appendix G) 

33.72 ( = C) 

S.D. = ((^(T - = 2.687 (°C) 

n=l 

where : 

S.D. = standard deviation 

Since logarithmic term in equation of "T^q' is too small 
when compared to standard deviation, this term can be 
neglected for uncertainty analysis (i.e., T^q = "^sat ^ ’ 



T = 33.72 (°C) 

wo 



6T 



wO 



2 X S.D. = 5.37 (°C) 



aT = 


T - T .= 33.72 - 

wo sat 


10.00 = 23.72 ( 


°c) 


6A T/A T 


= ((sicyAT)^ + (-5'^^t/AT)^)^ ' 2 




6 A T/ at 


= ((5.37/23.72)2 + (- 


0.5/23.72)2)1 ' 2 




6 A T/ at 


= 22.7 percent 






D. UNCERTAINTY IN HEAT FLUX 






q = 


(C^- 2Q,)/A^ 






3, = 


203.17 (W) 


6Q = 3.65 

H 

the uncertanity 


(W) 


Assuming 


the same proportion in 


for 




13.15 (W) 


6Q = 0.23 

1? 


(W) 


Ch - 


2Qp = 176.87 (W) 


r 
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6q/q = ( (6Q^/(Q^- ^ ^ 

(6A^/P^)^ )^ 

5q/q = ( (3.65/176.87)" + (0.46/176.87)" + 
(-0.0063 )" )' '" 

= 2.17 percent 

E. UNCERTANITY IN BOILING HEAT-TRANSFER COEFFICIENT 
h = q / A T 



6h / h = 


( 


(6q/q)^ ^ (-6AT/AT)" )^'" 


6h / h = 


( 


(0.0217)" + (-0.227)" )' '" 


- 


22 


. 8 percent 
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